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ABSTRACT 

The  backfire  bifilar  helical  antenna,  consisting  of  two  opposed  helical 
wires  fed  with  balanced  currents  at  one  end,  is  a  new  type  of  circularly 
polarized  antenna.  When  operated  above  thb  cutoff  frequency  of  the  principal 
mode  of  the  helical  waveguide,  the  bifilar  helix  produces  a  beam  directed 
alcng  the  structure  toward  the  feed  point.  The  term  "backfire"  is  used  to 
describe  this  direction  of  radiation  in  contrast  with  "endfire"  which  denotes 
radiation  away  from  the  feed  point. 

Radiation  patterns,  measured  for  a  wide.,  range  of  helix  parameters, 
show  maximum  directivity  slightly  above  the  cutoff  frequency.  The  pattern 

'  ”  “  f! 

broadens  with  frequency,  and,  for  pitch  angles  near  forty-five  degrees,  the 
beam  splits  and  scans  toward  the  broadside  direction. 

Near  field  measurements  show  the  current  decaying  rapidly  to  a  level 


-J 

dbout  twenty  decibels  below  the  input  level  at  a  rate  that  increases  with 
frequency.  Phase  measurements  in  the  near  field  show  that  in  the  feed  region 
the  direction  of  phase,  progression  is  toward  the  feed  point.  The  oppositely 
directed  phase  progression  and  direction  of  energy  flow  is  characteristic  of 
a  backward  wave.  The  direction  of  phase  progression  is  consistent  with  the 
backfire  direction  of  the  main  beam  observed  in  the  radiation  patterns  and 
the  increasing  rate  of  current  decay  is  consistent  with  the  broadening  of 
the  main  beam  with  increasing  frequency. 

A  theoretical  analysis  of  the  bifilar  helical  antenna  is  obtained.  It 


is  based  upon  the  semi-infinite  model  using  thin  wire  assumptions.  These, 

V .  .. 

so-called,  linearizing  assumptions  consist  of  replacing  the  current  distribu- 

«  ■  ■■■■  ....  " 

tion  on  the  surface  of  the  wire  with  a  line  current  onTithe  center  line  of  the 
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wire  and  of  satisfying  the  boundary  condition  along  one  line  on  the  conducting 
surface.  The  Fourier  transform  of  the  current  distribution  on  the  semi-infinite 
h^lix  is^ deduced  from > the  determine ntal  equation  of  the  bifilar  helical  wave- 
guide  by  a  Wiener-Hopf  techniqW/  The  relation  between  this  Fourier  transform 
and  the  radiation  pattern  of  the  backfire  bifilar  helical  antenna  is  shown. 

The  results  predict  the  patterns  of  the  experimental  study  and  s)iow  the  effect 
of  wire  size  on  antenna  performance. 
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1 .  INTRODUCTION 


1.1  Statement  of  the  problem 

The  purpose  of  this  study  is  to  provide  a  mathematical  model  of  the 

backfire  bifilar  helical  antenna  which  contains  the  essential1'  features  of 

■  '  -  n '  • 
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the  practical  antenna  and  which  is  amenable  to  solution.  This  model .has 


been  solved  for  the  far  field  radiation  pattern  of  the  antenna,  and  the 
results  are  compared  with  the  radiation  patterns  that  have  been  obtained' 
experimentally.  One  outstanding  characteristic  of  these  patterns  is  that 


e  principal  direction  of  radiation  is  in  the  backfire  direction.  The 


backfire  direction  of  a  surface  wave  antenna  is  along  the  structure  toward 


the  feed  point.  Radiation  in  this  direction  can  be  associated  with  a  phase 
progression  toward  the  feed  point,  while  the  group  velocity  must  be  away 
from  the  feed  point:.  This  is  characteristic  of  a  backward  wave  of  current 
on  the  structure.  The  term  "backfire-'  is  used  for  contrast  with  the  more 
usual  endfire  surface  wave  antenna  for  which  the  principal  direction  of  ra¬ 
diation  is  away  from  the  feed  point  of  the  structure. 

The  backfire  bifilar  helical  antenna,  shown  in  Figure  1,  ,1s  constructed 
of  two  helical  wires  wound  in  a  right  circular  cylinder  with  a  constant 
pitch.  The  corresponding  points  of  the  wires  are  located  at  the  ends  of  a 
diameter  of  the  cross  section  of  the  cylinder.  The  helical  wires  are  fed 
with  balanced  currents  from  a  source  on  the  axis  of'  the  cylinder.,  A  com¬ 
plete  description  of  the  antenna  includes:  thd\radius  b,  measured  from 
the  axis  of  the  cylinder  to  the  center  line  of  a  conductor;  the  pitch  dis¬ 
tance  p,  measured  along  a  generator  of  the  cylinder;  the  wire  radius  a, 
and  the  total  length  of  the  antenna.  .  ■if' 


1 .2  Method  of  Solution 

A  solution  of  ..the  backfire  bifilar  helical  antenna  .problem  is  obtained, 
first,  by  considering  |he  basic  helical  geometry  to  be  extended  to  infinity 
in  both  directions.  This  is  the  helical  wav.eguide'problem  which  already  has 
an  extensive  literature.  An  exact  account  of  the  boundary  conditions  in  this 
problem  involves  finding  a  current  distribution  over  the  surface  of  the 
conductor  which  makes  the  tangential  electric  intensity  vanish  everywhere 

on  the  surface  of  the  conductor.  An  approximation  used  by  Sensiper1  and 

2  . 

Kogan  is  to  require  that  the  tangential  electric  intensity  vanish  only 
along  some  line  alon^  the  surface.  This  approximation  is  good  for  rela¬ 
tively  thin  conductors.  Sensiper,  in  treating  the  tape  helix,  makes  a 
further  approximation  in  assuming  a  functional  dependence  for  the  distribution 
of  current  across  the  width  of  the  tape.  Kogan  makes  a  similar  approximation 
by  assuming  that ...  the  surface  current  distribution  may  be  replaced  by  a 
current  flowing  along  the  center  line  of  the  helical  conductor.  The  ap¬ 
proximations  involved  in  replacing  a  surface  current  by  a  line  current  and 
in  satisfying  the  boundary  conditions  along  only  one  line  of  the  surface 
have  been  termed  linearizing  approximations.  A  formulation  of  the  helical 
waveguide  problem,  similar  to  that  of  Kogan,  is  used  in  this  study. 

Using  the  above  approximations,  the  desired  current  distribution 
becomes  a  function  of  only  one  variable.  •  This  variable  may  be  either 
distance  along  the  helix  or  distance  along  the  helical  axis.  In  this  paper, 
the  helical  axis  is  taken  to  coincide  with  the  z-axis  of  a  cartesian 
coordinate  system,  and  z  is  taken  as  the  independent  variable  of  the  current 
distribution.  The  integral  equation  for  the  current  on  the  helical  waveguide 
now  has  the  form  of  a  convolution  integral.  This  happy  circumstance  related 


to  invariance  under  a  one  dimensional  Abelian  congruence  group,  suggests  a 
solution  by  Fourier  transformation.  The  determinantal  equation  for  the 
Fourier  spectrum  of  the  current  distribution  has  three  pairs  of  real  roots 
at  the  lower  frequencies.  As  the  frequency  is  increased,  two  pairs  of  roots 
tend  together  and  coalesce  at  the  critical  frequency,  f^.  This  frequency 
marks  the  lower  limit  of  the  range  of  frequencies  of  interest  in  this 
study.  The  critical  frequency  depends,  in  a  complicated  way,  upon  the 
radius  of  the  helix,  its  pitch,,  and  thtf ’conductor  radius,  and  no  simple 
expression  is  given  for  it. 

Tile  second  step  in  the  solution  of  the  backfire  bifilar  helical  antenna 
problem  is  to  obtain  a  solution  for  the  current  distribution  on  a  helical 
waveguide  extending  only  along  the  positive  z-axis  of  the  coordinate  system. 
This  is  accomplished  by  a  Wiener-Hopf  factorization  of  the  transformed 
solution  for  the  helical  waveguide.  The  boundary  conditions  on  the  factori¬ 
zation  are  such  that  the  current  is  identically  zero  for  all  negative  z  and 
approaches  a  finite  non-zero  limit  as  z  approaches  zero  from  the  right.  The 
discontinuity  in  the  current  at  the  origin  is  the  source  of  energy  in  the 
probiem|and  corresponds  to  a  pair  of  oscillating  point  charges  separated 
by  the  helix  diameter.  The  result  of  this  calculation  is  the  Fourier 
transform  of  the  current  distribution  on  a  semi-infinite  bifilar  helix  fed 
by  a  charge  dipole  in  the  plane  z  =  0, 

The  solution  obtained  in  this  way  is  limited  by  the  linearizing  approxi¬ 
mations  described  above.  The  remaining  approximations  in  the  theory 
developed  here  are  contained  in  equating  this  solution  to  the  radiation 
pattern  of  a  finite-length  bifilar  helical  antenna  fed  from,  a  source  on  its 
axis.  The  practical  antenna  differs1  from  the  mathematical  model  in  the 


il 
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vicinity  of  the  feed  point  by  the  addition  of  a  current  flowing  along  a 
diameter  of  the  cylindrical  cross  section.  This  'additional  current  element 
is  short  in  the  range  of  helixi<diam|€$r's  of  primary  •  interest  in  this  study 
and  will  not  contribute  noticeably  to  the  radiated  field.  Measurements 
of  the  current  on  several  models  of  the  bifilar  helix  above  the  critical 
frequency  show  that  the  current  decays  rapidly  with  distance  from  the 
feed  point.  At  a  level  about  20  decibels  below  the  input  value,  an  undamped 
wave  becomes  dominant.  Because  of  the  rapid  attenuation  of  current  on  the 

.•I 

bifilar  helix,  the  radiation  pattern  of  the  finite' structure  is  expected 
to  be  the  same  except  for  "end-fire"  radiation  caused  by  the  residual 
"free-mode"  wave.  It  is  well  known  that  a  finite  line  current  distribution 
and  its  far  field  radiation  pattern  are  Fourier  transform  pairs.  This  is 
extended  here  to  include  the  semi-infinite  structure,  and  it  is  shown  that 
in  this  case  also  the  radiation  pattern  is  simply  related  to  the  transform 
of  the  current  distribution  as  a  function  of  z  except  in  the  vicinity  of 
the  current  distribution  at  infinity. 

The  Fourier  transform  of  the  current  distribution  is  obtained  by 
numerical  techniques  using  an  automatic  digital  computer.  The  computed 
patterns  based  upon  this  Fourier  transform  show  good  agreement  with 
measured  patterns  presented  in  this  report.  This  confirms  the  validity 
of  the  approximations  used  in  constructing  the  mathematical  model  of  the 
antenna. 

1". 3  Review  of  Helix  Analysis 

The  study  of  electromagnetic  wave  propagation  on  helical  conductors 

3  1 

through  195,5  has  been  summarized  by  Sensiper  and  his  thesis  contains 
an  extensive  bibliography  of  the  literature.  He  obtains  an  exact  solution 


of  the  "tape"  helix  by  expressing  its  field  as  a  ..sum  of  sheath  helix  modes 

and  by,  requiring  the  tangential  electric. intensity  to  be  zero  on  the  tape. 

The  resulting  arithmetic  is  found  to  be  intractable,  and  he  quickly  obtains 

a  more  manageable  approximate  solution  by, requiring  only  that  the  tangential 

electric  intensity  at  the  center  of  the  tape  be  zero.  When  the  resulting 

2 

expression  is  compared  with  that  obtained  by  Kogan  from  a  potential 

integral  formulation,  they  are  found  to  differ  principally  in  the  form  of 

the  convergence  factor  introduced  by  the  approximation.  Kogan  treats  the 

helical  wire  model  with  the  null  field  boundary  condition  satisfied  only  at 

the  points  of  tangency  between  the  helical  conductor  and  a  cylinder  with 

radius  equal  to  the  outer  radius  to  the  helical  wire. 

The  work  of  these  authors  is  directed  toward  the  evaluation  of  the 

'^free-mode"  propagation  constants  for  undamped  traveling  waves  of  current 

on  the  helical  conductors.  Although  Sensiper1  does  devote  some  space  to 

the  source  problem,  his  formulation  of  this,  problem  is  for  the  infinite  helix. 

It  is  used  to  interpret  the  "free-mode"  solutions  of  the  source  free  problem. 

Many  authors  have  contributed  to  the  literature  on  the  helical  waveguide. 

1,  3, 

Since  thbir  work  has  been  reviewed  by  Sensiper  ’  ’  the  bibliography  is  not 

reproduced  here , 

l.fl  Helical  Antennas 

The  helical  antenna  which  most  nearly  approaches  the  backfire  bifilar 

helical  antenna  in  performance  is  the  helical  beam  antenna  introduced  by 
'4 

Kraus  in  1947,  Kraus's  antenna  is  a  monofilar  helical  wire  fed  at  one 

ii 

end  against  a  ground  screen.  The  properties  of  this  antenna  have  been 

discovered  by  experimental  techniques.  The  current  distribution  on  helical 

"  "  5 

structures  of  this  type  was  studied  by  Marsh  .  The  analysis  of  a  helical 
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antenna  of  this  type  suffers  the  same  difficulties  as  the  analysis  of 

the  usual  end-fire ^surface  wave  antenna.  For  both  of  these  the  feed 

■v\ 

region  and  the  terminating  region  are  important  in  determining  their  radiation 
characteristics.  In  contrasty  the  characteristics  of  the  backfire  bifilar  . 

helical  antenna  are  controlled  by  the  input  region. 

,  "  6 

A  bifilar  version  of  the  Kraus  helix«nas  been  reported  by  Holtum  . 

This  antenna  is  constructed  of  two  coaxial  helical  wires  on  the  diameters  of 
the  supporting  cylinder.  Each  conductor  is  fed  against  a  ground  screen  with 
the  exciting  currents  in  phase ^opposition.  This  differs  from  the  backfire 
bifilar  helical  antenna  in  which  the  conductors  are  fed  at  one  end,  one 
against  the  other,  without  the  presence  of  a  ground  screen.  It  is  the 
absence  of  the  ground  screen  that  distinguishes  the  backfire  helix  from 
the  Kraus-type  helix.  This  fact  is  essential  to  the  performance  and 
analysis  of  the  backfire  helical  antenna..  On  the  other  hand,  the  number 
of  helical  conductors  is  not  essential  to  the  backfire  characteristic 
of  the  antenna.  A  backfire  monofi^ar  helix  is  shown  to  have  substantially 
the  same  radiation  characteristics  "as  the  backfire  bifilar  helix.  The 
monofilar  helix,  however,  is  more  difficult  to  feed  in  'the  backfire  mode. 

Although  both  the  Kraus  helix  and  the  backfire  bifilar  helix  are 

'S 

w 

circularly  polarized,  several  differences  in  the  performance  of  these 
antennas  must  be  noted.  The  beam  width  of  the  Kraus  helix  decreases  with 
frequency  while  the  backfire  helix  beam  width  increases  with  frequency. 

The  gain  of  the  Kraus  helix  increases  with  length  while  the  gain  of  the 
backfire  helix  is  independent,  of  length  provided  the  length  is  large 

-  ..  .  ..  "  ..  ..  H  " 

enough.  Finally,  the  Kraus  helix  is  an  end-fire  antenna  in  contrast  to 
the  backfire  helix  which  radiates  along  the  helical  axis  toward  the  feed 
point  r  ... 


w 


«  •  8  . 

1.5  Organization 

In  the  present  section,  the  backfire  bifi lar» hell  cal  antenna  is  described 
and  compared  with  other  helical  antennas  of  similar  character.  A  description 
of  the  method  used  to  obtain  a  solution  for  the  radiation  pattern  of  the 
antenna  is  given,  and  a  brief  review  of  the  literature  on  ^the  helical  wave¬ 
guide  is  included.  In  Chapter  2  the  relation  between^ 'tj£#jf radiation  pattern 

..  ..  v  rf " 

of  a  helical  antenna  and  the  Fourier  spectrum  of  its  current  distribution 
is  discussed.  This  motivates  the  work  of  the  three  following  chapters  in 
which  the  Fourier  spectrum  of  the  current  distribution  is  obtained. 

Chapter  3  describes  the  approximate  dete rminantal  equation  for  the 
helical  waveguide  as  it  is  used  in  this  study.  This  is  derived  from  the 
potential  integral  with  a  linear  current  approximation.  The  equation 
describing  the  helical  waveguide  is  factorized  in  Chapter  4  by  the  method 
of  Wiener-Hopf,  and  the  evaluation  of  the  resulting  expression  by  numerical 
techniques  is  described  in  Chapter  5,  The  results  of  the  numerical  computation 
are  presented  in  Chapter  6  along  with  the  experimental  results.  This  Chapter 
also  indicates  the  connection  between  this  study  and  the  log-spiral  and  other 
frequency- independent  antennas.  Chapter  7  summarizes  and  concludes  the  work. 


2.  THE  HELICAL  ANTENNA  RADIATION  PATTERN 


2.1  Formulation  of.  the  Vector  Potential 

Since  the  purpose  of  this  study  is  to  develop  a  theory  that  is  capable 
of  predicting  the  radiation  patterns  of  the  backfire  bifilar  helical  antenna, 
and  since  the  geometry  of  the  problem  lends  itself  to  a  solution  by  Fourier 
transformation^  a  relation  between  the  radiation  pattern  and  the  Fourier 
spectrum  of  the  current  distribution  for  the  bifilar  helix  must  be  obtained. 
This  section  is  devoted  to  establishing  that  relation. 

The  radiation  pattern  of  an  antenna  will  be  taken  to  mean  the  angular 
distribution  of  electric  intensity  in  spherical  coordinates  at  a  large  dis¬ 
tance  from  the  antenna.  When  only  those  teriiis  that  vary  inversely  with  dis¬ 
tance  are  retained^  the  electric  intensity  is  related  to  the  magnetic  vector 
potential  by 

£xr  =  -tkC  Axr  (1) 


where  r  is  a  unit  radial  vector 


k  =  u  is  the  propagation  constant  of  free  space 

is  the  intrinsic  impedance  of  free  space 
Thus  the  radiation  pattern  is  related  to  the  angular  distribution  of  the  vector 
potential,  and  the  vector  potential  is  deduced  from  the  current  by  the  well 
known  integral  formula 


A(r)  =  r  G(r,  r  )  I (r  )  dv 


where 


10 


is  the  Green's  function  in  an  unbounded  homogeneous  isotropic  three  dimensional 
space," 

..  \\  •>..  . 

r  is  the  radius  vector  to  a  point  of  observation  and 

—  V/ 

r  is  the  radius  vector  to  a  source  point. 

By  the  method  of  induced  sources,  we  may  replace  the  conducting  boundaries 
at  the  surface  of  the  helical  wires  by  a  surface*  current  distribution  existing  » 
in  homogeneous  space.  If  this  current  distribution  equals  the  surface  current 
on  the  helical  conductors,  the  fields  external  to  the  conductors  are  unchanged 
by  removing  the  conductors  from  the  space.  If  the  conductors  are  sufficiently 
thin,  the  field  produced  by  the  surface  current  distribution  will  not  differ  much 
from  that  produced  by  a. current  distribution  on  the  center  line  of  the  helical 
wire.  This  approximation  limits  the  analysis  to  thin  wire  helices,  but  it 
has  the  advantage  of  converting  the  volume  integral  in  Equation  (2)  to  a  line 


integral. 


The  bifiiar  helix  then  can  be  considered  as  two  helical  lines  of  current. 

/ 

A  point  on  line  one  with  axial  coordinate  z  will  have  cartesian  coordinates 

P  -  (b  cos  Tz/",  b  sin  Tz/,  z/) 

where 

t  =  2ir/p 


p  =i',  pitch  of  helix 
b  =  radius  to  centerline, 

and  the  unit  tangent  to  the  helical  wire  at  that  point  is  given  by 

A'  4  "  i,  /  ^  I  /  A  ■  ■  I 

-x  cos  4^  sin  Tz  -fy  cos  4*  cos  Tz  +z  sin  41 

where  41  is  the  pitch  angle  of  the  helix  given  by 

■  i'v"  ,  27Tb  ,  - 

-,co.t  ^  - ■■■■  Tb  » 

P 


XI 


At  the  same  value  of  the  axiai  coordinate  z*a  point  on  line  two  will  have 
cartesian  coordinates 


P  =  (-b  cos  Tz*,  -b  sin  rz' }  z' ) 


with  unit  tangent 


"  —  x  cos  ^  sin  Tz  -  y  cos  41  cis  Tz/  +z  sin  41  .. 

The  distance  between  a  remote  point  (r,  0,  <f>)  and  a  point  on  line  1  is 


given  by 


r 

=  r^  =  (r  sin 


v  /  2 

in  0  cos  tf*-*3  cos  tz  ) 


4  (r  sin  0  sin  4>-b  sin  Tz  )  +  (r  cos  6 


->2] 


>  -  P- 


2br  sin  6  cos  (Tz  -<J>)  4 


u2  2”] 
b  +z 


Similarly,  the  distance  between  the  remote  point  and  a  point  on  wire  2  is 


given  by 


L2  2 

r  +2br  sin  0  cos  (Tz-cj>)  fb  +z 


The  currents  on  the  helical  line  can  be  considered  to  be  a  function  of 
the  axial  variable  z/  alone.  If  the  currents  in  the  two  wires  are  equal  and 
oppositely  directed  for  the  same  value  of  z Equation  (2)  may  be  written 


Ax(r)  =  -  cos  4*  J  +  '4i Ff  )  Sln  TZ'  1 

r)  =  cos  *  ‘  C°S  TZ'  1  (z/) 

')  =Sln'  j  ~  ^rr~)  1  (z'5  dz' 


I  (,z/ )  dz 


T  /e-Jkri  e-J*r2 

AxCr)  =  cos  +  j  -i;  ----- 


I  (z')  dz' 


2  .-2  Transformation  by  Parseval 1  s  Theorem 


In  Equation  (3),  the  vector  potential  at  the  point  (r,  Q,  <p)  can  be  related 


to  the  Fourier  spectrum  of  the  current  distribution  by  Parseval's  theorem. 
Parseval's  theorem  states:  Given  two  functions,  F(z)  and  I(z),  of  the 

r>J 

spatial  variable  z  with  Fourier  transforms,  F(P)  and  I(P),  in  the  transform 
variable  P,  " 


F*  (z) 


I.(z)  dz  =  27T 


f(P)..i(P)  dp 


where  I(P)  is  related  to  I(z)  by 

Tcp)  "  ^  L uz)  gjPz  dz 

Koo 


This  introduces  the  Fourier  transform  of  the  current  distribution  into  the 
calculations  for  the  vector  potential.  In  this  paper  the  tilde  (~)  is  used 
to  ..denote  the  Fourier  transform  of  a  function  of  z,  and  the  asterisk  (*)  is 
used  to  denote  the  complex  conjugate  of  a  function.  It  is  now  necessary  to 
find  the  Fourier  transform  of  the  remaining  factor  in  the  integrand  and  to 
evaluate  the  transformed  integral  at  a  point  remote  from  the  helix.  The  complex 


conjugate  of  this  factor  in  the  first  integral  of  Equation  (3)  can  be  written 

Jkr, 


F*(zj 


-  /cJ 


2  eJkr2 
+ 


sin  tz  (gj+gg) 


ejTZ-e-jTz 

8ffj 


where  g  and  g  are  the  functions  discussed  in  Appendix  A  i 

“  2-  2  2  1/2 
exp  (jk[A  +2AB  cos  (6-frT z)  +B  +(z-d)  ]  ) 


S-, 


2  _  2 

[A  +2AB  cos  (tjji’Tz)  +B  +(z-d) 


2,1/2 


where 


A  =  r  sin  0 
B  =  b 

d  =  r  cos  0  £ 

Using  Equation  (A-5)  and  the  shifting  theorem  for  Fourier  transforms  we  obtain 


F*(B)=  777  E 

16TT 

n  even 


J  (b[k2~(|3-(n-l)T)2]1/2)  H(1)(r  sin  e[k2-(B-(n-l)T)2]1/2) 

1  11  5  .  n 

j(B-(n-l)T)  r  cos  6  +.jn<j> 

.  e  e 

/-J  (,b[kZ-((B-Cn+l)T)2]1/2)  H" Cr  sin.  0[k2-(B-(n+l)T)2]1/2) 


j(B-(n+l)T)  r  cos  6  jn4> 


In  this  way,  using  Parseval's  theorem,  the  first  equation  in  (3)  becomes 

,2  ,  _ 2,1/2.„(1),  or,  2 


00 

.  ,  X  +COS  ^  ~  , 

Ax(r)T“E  -00  Wb] 

n  even  1 


2nl/2x„(D  ~ j  -  fir tn  /«  2nl/2x 

n 

j(B-(n-l)T)  r  cos  9  jn<j> 


(j  (b[k'i-(B-(n-l)T)'5]1^)H^y(r  sin  0{k  -(B-<n-l)T)*] 
n  n  .. 


-J  (b[k2-(B-(n+l)T)2]1/2H(1)(r  sin  9[k2-(B-(n+l)T)2];l'/5 
\  n  n 


ej(B-(n+l)T)  r  cos  sejn<p 


2.3  The  Far-Field  Radiation  Pattern 


The  asymptotic  estimation  of  this  integral  follows  that  given  in  Appendix  B. 
In  this  case  the  change  of  variables  (Equation  B-4)  is 

6  -(n-1)  T  =  k  cos  a 


or 


Thus  we  have 


B  -(n+1)  t  =  k  cos  a 


A  (r)~+j7T  cos  4*  E  G  J  (kb  sin  Q)[I(k  cos  0+(n+l)T)-I(k  cos  0+) (n+l)T) ]e' 
x  „  o  n 

n  even 


Jn't’ 


where 


-jkr 

G  =  e  ° 
o  4TTr 

o 


In  this  way  Equation  (3)  becomes  // 

A  ("ri  ~  +j7T  cos  41  £  G  J  (p)[I(k  cos  04:(n-l)T)  -l(k  cos  9+(n+l)t)  ]e^n<^ 


(4) 


A  (r^  ~  JTcos  4*  £  G.J  (P)[l(k  cos  0+(n-l)T)  +I(k  cos  9+(n+l)T)  ]e^n^ 
J  n  even  °  n 

A  (r)  ~  2 ft  sin  4  L  G  J  (P)  I(k  cos  04-n'Oe^1^ 

z  -  •  ,•  ,  ,  o  n  ,  .  .. 

n  odd  ••  - 


if 


II 
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where  „  ,,  ■ 

P  =  kb  sin  0 

~  4>\ 

These  series  are  rapidly  convergent  if  I(k  cos  9+nT)  is  bounded  on  n  because  of 

the  properties  of  J  (P) .  The  ratio  of  the  first  two  of  these  to  J  (P)  is  shown 

n  o  . 

in  Figure  2.  In  addition,  it  will  be  shown  later  that  for  the  backfire  bifilar 

helix  |l  (k  cos  0+nT)  j  also  decreases  rapidly  with  n.  Thus  the  radiation  field 
n 

of  the  backfire  bifilar  helix  may  be  represented  by  the  first  terms  of  Equation  (4) 
with  little  error  giving 

A  (r)  ~  J7T  cos  f  G  J  (P)[I(k  cos  0-T)-i(k  cos  9+t)  ] 


A  (r)  ~  TTcos  f-G  J  (P)[I(k  cos  6  +T)  +  I(k  cos  0-T)  ] 
y  o  o 


(5) 


~  27T  sin  t  G  J  (p)[I(k  cos  8+T)e^+  I(k  cos  0-T)e  JT] 
z  o  1 


From  Equation  (5)  it  is  obvious  that  the  backfire  bifilar  helix  radiation 
pattern  is  the  sum  of  right-handed  and  left-handed  circularly-polarized  components. 
The  ratio  of  these  components  on  the  axis  of  the  antenna  is  given  by 


r  = 


I(l-k) 


I(il-k) 

The  axial  ratio  of  the  polarization  ellipse  is  deduced  from  this  by  the 
equation 


AR  =■ 


1+r 

1-r 


(6) 


Computed  values  of  the  axial  ratio  are  described  in  Chapter  6  of  this  report. 

The  relation  between  the  radiation  pattern  of  the  backfire  bifilar  helix 
and  the  Fourier  spectrum  of  its  current  distribution  is  established.  The 
spectrum  of  the  current  distribution  on  the  semi-infinite  helix  will  be  deduced 
from  the  determinantal  equation  for  the  infinite  helix.  The  determinantal  equation 
for  the  infinite  bifilar  helix  is  derived  in  the  next  chapter. 


CIRCUMFERENCE  IN  WAVELENGTH  P 

Figure  2.  Ratio  of  Higher  order  terms  to  the  First  Term 
in  the  Series  Expansion  for  the  Pattern  of  the 
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3.  THE  BIFILAR  HELIX  DETERMINANTAL  EQUATION 


(7)| 


3.1  The  Complete  Circuit  Equation 

•  ••  //  . 

The  electromagnetic  field  in  an  infinite  homogeneous  isotropic  medium  can 

be  obtained  from  the  magnetic- vector  potential  A  and  the  scalar  potential  V 

/ 

using  the  relations 

H  =  VxA  -  |f 

I!  , 

£  =  „-jk£  A-W 

-k. 

If  the  divergence  of  A  is  chosen  to  be 

V-"a  =  -jk/£  V 

these  potentials  satisfy  the  scalar  and  vector  wave  equations, 

fl^+k^V  =  -~ 

=  -  ? 

The  solution  of  these  equaltions  can  be  written  in  terms  of  the  Green's 
function  of  the  medium 


exp[-jk|r-^|  ] 
G(rj  r  )  =  47T  - 


V(r)  =  z  ~  )P(*  )  dv 

A(r)  =JJJa(r,  ?0  I(r')  dv 

Using  the  equation  of  continuity 

V  -I  =  rjw p 

in  Equation  (8)  with  Equation  (7),  there  results 

V  =  %/k I  I  h(r,  P)  V-I  (P)  dv 

W  =  -E  -jk&  I  llG(r,  rO  I(r)  dv 
which  forms  the  basis  of  the  complete  circuit  equation. 


(8) 


(9) 
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The  circuit  equation  for  the  hlf liar  helix  is  obtained  by  replacing 

the  conducting  boundaries  of  thehelical  wires  by  the  current  distribution 

..  ((..  '  .. 

oh  the  surface  required  to  produce  dero  tangential  electric  intensity  there. 

•  ■"  '  i) 

The  equation  is  simplified  by  “Linearizing  the  boundary  condition.  This  is 
based  on  two  thin  wire  approximations;  first,  it  is  assumed  that  the  current 
distributiop  on  the  surface  of  the  wires  can  be  replaced  with  a  line  current 
distribution  at  the  center  of  the  wire,  and,  second,  it  is  assumed  that,  if 
the  electric  field  is  zero  along  only  one  line  of  the  surface,  the  resulting 
solution  will  be  a  reasonable  approximation  to  the  exact  solution.  These 
■assumptions  are  good  if  the  wires  are  sufficiently  thin. 

In  this  analysis  the  current  will  be  assumed  to  exist  on  the  line 

•\\ 

defined  in  cartesian  coordinates  by  ..  • 

p^  =  (b  cos  Tz/[  b  sin  Tz/,  z/) 
p  =  (-b  cos  Tz/,  -b  sin  Tz/,  z/) 
and  the  null  line  will  be  taken  as  defined  by 

q^  =  (b/  cos  Tz  }  b  'gin  tz,  z) 

q2  =  C°S  TZ>  Tz;  z ) 

where 

b  *  -  b-a  "  • 

II 

b  =  distance  from  helix  axis  to  wire  centerline 
a  =  wire  radius 

Thus,  the  currents  and  potentials  can  be  considered  a  function  of  the  z 
coordinate  alone.  Because  of . the  constant  pitch  angle  of.,  the  helix, 

=  cot  1  (27Tb/p) 

an  elementary  displacement  along,  the  axis  is  related"  to  one  along  the  wire  by 

dz  =  sin  4*  dp  .. 
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Thus  on  wire  one 


V-I(r')  =  sin  + 


dl(z') 


W(r')  =  sin  ^  u  (z) 

1  dz 

I.(r ' )  =  u^z')  I^z') 


where  =  -x  cos  ^  sin  Tz+y  cos  cos  Tz+z  sin  ^ 

In  these  equations  it  is  assumed  that  the  pitch  a.ngle  of  the  null  line  *lj/is 
equal  to  the  pitch  angle  4*  of  the  centerline  of  the  wire.  The  correct  result  is 


where 


cot  L|j/  =  —  cot  +  =  (1-&)  cot  + 

D 


5  =  a/b 


The  approximation  OP  =  +)  is  consistent,  with  the  thin-wire  assumption  used 
throughout  this  study. 

3.2  The  Determinants!  Equation  for  the  Monofilar  Helix 

Equation  (9)  will  be  applied  first  to  a  single  wire  helix  so  that  the  , 

modification  of  the  determinants!  required  for  the  balanced  bifilar  helix 

will  be  made  evident.  The  current  is  confined  to  the  line  and  the  potential 

is  evaluated  on  the  line  q^  Since  the  electric  field  strength  is  assumed  zero 

along  this  line,  Equation  (9)  becomes 

do  d_ 

v(z)  =  j  |  J  Gn(z>  z  )  sln  ^  dz~ ~  dz' 

-00 

. 00 

dv  r 

sin  +  =  -jk4  J  G^^(z,  z')  u^CzJ-^fz')  I^z')  dz '  (10) 

-oo 

If  the  first  equation  is  differentiated  with  respect  to  z,  multiplied  by  sin 
and  equated  to  the  second  equation;  the  result  is 


[sin  4* 


dG  (z,  z  ) 


dl<z0  ,  ,2  .a 


+  k  (u^z^u^z')  G1]L(z,  z')  I(z ')  ]  dz  (11) 


is 


G  (B)  E  J  (b[k2-(B-nT)2]1/2)  H<2)(b[k2-(B-nT)2]1/2) 

11.  o  (/  ...  _  n  ■  n 

#  n=-oo 


x 


I„ <b [ (B-nT) 2-k2 ] 1/2>  Kn (b[ (B-nT) 2-k2 ] 1/2) 


=  77  E  n 

u  4  7]^  n=-tt 

For  convenience  we  will  define  a  new  function 

B  (B)  a  I  (b/[(B-nT)2-k2]1/2)K  (b[(3-nT)2-k2]1/2) 
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The  function,  0^(6)  now  has  the  simpler  representation 

'  Lj‘w- 


The  Fourier  transform'of  Z(z)  is  given  by 

k2cot2^  r~ 


Zn(l3)  =  <B2-k2)  Gu(B) 


[Gu  (B+T)+G11(B-T)] 


This  function  has  infinitely  many  branch  points  located  at 


B  =  nT±k 


The  Fourier  transform  of  Equation  (12)  is 


n  =  0,  -X,  -2,  •••■■ 


ZU(B)  I X(B)  =  0 


(13) 


The  spectrum  of  the  current  on  the  infinite  monofilar  helical  waveguide  given 

by  Equation  (13)  must  be  a  distribution  of  point  support  at  the  roots  of 

Z(B)  and,  therefore,  is  a  polynomial  of  delta  distributions  there.  However,  by 

7 

analytic  continuation  of  the  distribution  into  the  complex  B-plane,  it  can  be 
shown  that  the,  order  of  the  distribution  must  be  one  less  than  the  order  of  the 

/V> 

root.  Since  the  roots  of  Z(B)"  are  of  first  order,  no  derivatives  of  the  delta 
appear.  The  current  distribution,  therefore,  is  a  sum  of  traveling  waves  with 

rv 

propagation  constants  given  by  the  roots  of  Z(B). 

<v  2 

In  investigating  the  roots  of  Z(B),  Kogan  rearranges  terms  in  the  form 


i.  n 


"  ..  "  ..  *  '  ■  ..  .  /  '  ' 

ii 

oo 

S  [B  (p+T)  +  B  (p-T)] 

n--°o _ _ _ ■■ 

00 

2  2  B  .  (6) 

n  ’ 

n — °° 

2 

The  function  on  the  left  can.be  called  F^(p)  and,  following  Kogan  ,  is  sketched 
in  Figures  3  and  4,  It  can  be  seen  that  F^fp)  aPProximately  equal  to  one  in 
the  range 

|  <  I  |  -  n  I  <  1  -  £  n  =  0,  +1,  +2,  .  .  .  ' 

“4* 

except  near  p  =  r  —  k  where  the  terms  B^(p-T)  in  the  numerator  have  a  logarithmic 

singularity  contributed  by  the  modified  Hankel  function  of  zero  order  or  near 

P  =  k  where  BQ(p)  in  the  denominator  is . logarithmically  singular.  In  this 

range  the  root  of  Z  (p)  is  given  approximately  by 

P2  2  i 

•E—  =  1  +  cot  41 

.  k2 
or 

p  =  k/sin  + 

The  roots  of  Z(p)  are  at  the  intersection  of  F  (p)  and  the  parabola 

A 

\ 

For  k  small  it  can  be  seen  that  there  are  three  roots:  one  near  p  =  k/sin  41, 
one  near  p  -  1-k,  and  one  near  p  -  1+k.  Since  Z(p)  is  an  even  function  of  p, 
there  are  also  three  corresponding  roots  on  the  negative  p-axis  . 

The  functional  dependence  of ' the  roots  of  Z(p)  on  the  frequency,  k,  is 
conveniently  displayed  on  the  Brillouin  diagram,  also  called  the  p-k  diagram, 
shown  in  Figure  5.  In  this  figure,  due  to  Sensiper1,  the  variables  are 
normalized  with  respect  to  t  (T_:2F/p).  The  deviation  of  the  curve  from  the 
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asymptote, 


..  p  =  k/sin  4 

near  p  =  0  can  be  interpreted  as  coupling  between  the  current  wave  and  a 

linearly  polarized  plane  wave,  and  the  deviation  near  p  =  .T+k  can  be  interpreted 

as  coupling  with  circularly  polarized  plane  waves.  This  coupling  is  indicated 

by  the  logarithmic  singularities  and  depends  upon  the  wire  thickness  6(6=a/b) 

of  F  (p)  shown  in  Figure  3. 

2  « 

3.3  The  Determinantal  Equation  for  the  Bifilar  Helix 

When  the  second  conductor  is  added  to. the  problem,  Equation  (12)  becomes 

oo  "  "  .. 

f  [ Z^1  (z-z 1  )  I^z’)  +  Z12(z-z*)  I2(z’)]  dz*  (14) 


where 


Z12(z_z,)  ;■  G12(Z"Z,)  h  ■  k2ci2(z"z,)  G12(Z-Z,) 


In  this  equation 


12 


22  2  1/2 
,  , ,  exp(-jk[b  ib'  f2bb'  cos  t(z-z' )+ (z-z1 )  ]  ) 

(Z-Z  )  .....  _  -  o“T79 - 

471  [b  tb'  +2bb'  cos  t(z-z’)» (z-z')  r/z 


relates  the  potential  at  wire  one  to  the  current  on  wire  two,  and 

C  (z-z')  -  u  (z)  •  u  (z')  ~  -  cos2  4  cos  T(z-z’)  +  sin2  4* 

Lei  1  4 

is  the  cosine  of  the  angle  between  the  line  element  at  the  source  point  and 

V 

the  line  element  at  the  observation  point. 

The  fourier  transform  of  G  (z)  is 


G12(P)  £  (-!)11  B  (?) 

471  n--.:-°c 


and  that  of  Z12(z.)  ls 


zncp)  =  cp2-  k2:  o12(p)  +  ---y  ?  [312(p+T)  +  g12(p-t)] 
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0^(8)  - 
11  27 T  -00 


00 

a-.  J 


-jk[b2+bf2-2bb/(/< cos  Tz+z2  ] 1//2 
47T[b2+b/2-2bb  cos  TZf.z2]1/'2 


ejBzdz 


'  \i  00 

-V  B  (8)  +  L  (B  (B)+B  (-B)) 

°  .2  „  '  n  n 

417  4it  n=l 


(18) 


where 


B  (B)  =  I  (b'vn)  Kn  (b  yn) 
n  ‘  n 


V  2  =  (3-nT)2  -k2 
n 


v# 

/.> 


When  k<T  all  the  terms  in  the  series  Equation  (18),®<are  positive  except  at 

most  one  term.  It  follows  that,  it  the  series  representation  of.G  (B)  converges 

then  G  (B)  converges  absolutely. 

1.2 


A  result  of  Watson 


8 


oo 


I  (x)  K  (x)  ~  I  — § — —  j2  (z)  dz 

n  n  J  ^+x2  n 

o 


shows  that  I^fx)  K^Cx)  ls  positive  monotone  decreasing.  Since 


Urn  I  (x)|K  (x)  =  — 
x->o  n  ▼  n  2n 


it  follows  that 


I  (x)  K  (x), 
n  n  ” 


<  1_ 
2n 


for  all  real  x.  The  recurrence  formula  for  the  modified  Hankel  function 

dxflK  (x) 

—  "  =  -x11  K  (x) 

dx  n-1. 

tells  us  that  xnKn(x)  is  also  positive  decreasing  for  positive  x  and,  if 

x,  >x  .  it  follows  that 
1  2* 


Thus 


x  K  (x:  )  <  x  *K  (x  ) 
I  n  I,  '2,  n  2 


*2  n  V  '  1  xo  n 

I  (x  )  K  (x  )  <  (-*)  I  (x  >  K  (x  )  <  —  (~S) 

n  2  nix  n  2  n  2  2n  x, 

...  A  5  ■  '■  .  1 


exp  <[n2+(b  Y  )2]1/2-n  sinlT1-^ 


Kn(bV  >■ 

Hi 


bv„ 


[ n2+(by  )2] V4 


The  term 


[n2+(bvn)2]1/2  =  [  ti2yi2(bT)2-2n  ®  (bT)2+  f(|)2.(^)2](bT)2  ]1/2 
can  be  approximated  tor  large  n  by 


[  „2+(bv>2]1/a  S  „  wWV  -  2  ££i2 


fcot2^ 


n  3  cos  41 

sin  4*  t  sin  41 


where  the  relation  .. 


bt  =  cot  41 


has  been  used.  The  term 


n  .2,1/2, 


n  sinh'1  -  n  in  <-L  +  [l+<-S*)a]V*) 


can  be  written 


-1  n  /n+[n2cot24'-2  |  n  cot2^  -  (£)2cotV+n2  J1/2 

n  sinh  —  =  n  In  - - - -  T _ _____ 

Yn  \  cot  t(|  -n)2-(£)2]1/2 


and  approximated 


n  sinh  1  — —  ~  n  in)--8i-  - _ I  sin  41 

I  n  cot  41  (1-  —  ) 

V  »  nT 
i,' 

n  sinh"1  ~  n  in  [  i_  !L  cob2*  ,  r.  B 


n  B  cos2^1 


bV 


cos 


T-  1  nT  uiTHT  ]  [1+  i?l> 


r: 


~  r*  in  [i-  3-SAP  *  1) 

.  Gos.  y- /" 1  ■■ 


Thus  the  modified  Handel  function  can  be  approximated  by 


l 


K  (by  ) 
n  n 


exp^  " 


B  cos  4*  -I  n  (> 

f  T  sin  ^  /  l+;sin  ^  ,  ,  ,  8  .  ,i,  \ 

—  (  -c^r>  ^  1+nT Bln  ^ 


but  since, 


,  ,  /,  8  sin  4i  \  8  sin  41 

lim  (It — )  =  exp 

n—^oo 


we  have 


:  8 


ir  sin  ^  , 1+sin  4\  2  n~T 


K  (bY )  ~  IT. 1^5 

n  n  V  2n  cos  111 


eXpC- 


Similarly,  using 


b  *T  =  cot  H'  =  (1-6)  cot  41 


we  have 


.  x  /  sin  ^  .cos  V*  k  r  n-8/T 

vb  V  'V  e^p  c 

\ 

Thus  the  Bessel  function  product  in  “Equation  (18)  can  be  approximated  for  n 
.large  by 

„  yn\  /n>  </ S  t «  +  Sill  >K  fcOS  ^(l+Sln  "  ,  3,  ,  1  1  x 

n  n  2n  Jcos  t  (1+sin  Y  J  T  sin  Y  sin  Y 

The  rate  of  convergence  of  the  scries  representation  of  g.^(3)  can  be  increased 

by  using  Hammer's  transformation10.  Since  the  terms  B/(S)  can  be  summed,  the 
U  n 

result  Is 


(8)  =  i  |-Vsln  ^  8ln  ^  cosh  ? 

4TT 

cos  4*  (1+sin  40 


11 
in  [l  - 


cos  4*  (1+sin  41  ) 
oo 


exp  ^  sin  4J  sin  4*' 


-)] 


+B  (B)  +  £  [B  (8)-B  (8)+B  (-B)-B  (-8)1 
o  ,  n  n  n  n 

n=l 


(19) 


Equation  (19)  is  in  the  form  used  in  the  numerical  calculations  of  this 
study..  The  expressions  in  the  first  term  are  somewhat  complicated,  and  it- 
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is  interesting  to  apply  further  approximations  based  on  the  assumption  of 
small  relative  wire  size  6.  The  algebra  involved  is  straightforward  but 
tedious.  The  result  for  the  first' term  is 

6 

-sin  +  In  — — tt  (20) 

sin  y 

This  term  controls  the  "sharpness"  of  the  corners  of  the  Brillouin  diagram. 

It  determines,  in  effect,  the  degree  of  coupling  between  the  current  wave  on 
the  helical,  wire  and  the  axial  plane  waves  of  linear  and  circular  polari¬ 
zation.  As  this  term  beoomes  larger,  the  corners  become  sharper.  When  it  is 
small,  the  corners  will  be  well  rounded,  and  the  coupling  can  be  considered 
large.  The  magnitude  of  this  factor  is  controlled  by  the  relative  wire  size, 

6,  and  the  sine  of  the  pitch  angle.  As  the  wire  size  increases  and  also  as 
the  pitch  angle  decreases,  coupling  is  increased,  and  the  Brillouin  diagram 
will  depart  considerably  from  the  asymptote 

8  =  k/sin 

On  the  other  hand,  when  the  wire  is  thin  and  the  pitch  angle  is  large,  coupling 
is  small  and  the  asymptotic  form  of  the  Brillouin  diagram  well  approximates  the 
current  distribution. 

It  proves  to  be  quite  difficult  to  obtain  an  asymptdtic  estimation  of 
Qn(B)  for  large  B  from  the  series  representation.  It  can,  however,  be  obtained 
rather  easily  from  the  integral  representation  in  Equation  (18).  The  integrand 
has  branch  points  at  the  roots  of  the  denominator 

2  12  r  2 

b  +b  -2bb  cos  Tz-t-z  -  0 

This  equation  has  no  real  roots.  If  z  =  x-f-jy,  the  real  and  imaginary  parts 

'  *  -  '  --  '  ,  )  ■“ 

of  this  equation  become 
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the  y  axis.  The  contribution  of  the  branch  point  at  z  =  j  1)/ t  is  evaluated 
by  the  method  of  van  der  Waerden11'.  According  to  this  method 

■  •  I/O 

S  =  j  ( jz+n/T) 

is  chosen  as  a  local  uniformizer  and  a  Laurent  expansion  of  the  integrand  about 
S  -s=  .0  is  obtained.  ,  " 

expC-jktb^+b7  2-2bb/  cos  tz4.z2]1/2)  a_l 

i t - - - — - - — i - —  =  -  +  a  +  a  S  +  .  .  . 

r,2  ,  2  /  _  2,1/2  S  o  1 

[b  +b  -2bb  cos  Tz+z  J 


where 


-J 


_1  V -2(Tbb/  sink  2/T[  (l-6)  cot^  sink  ??+??] 

The  integral  Of  the  first  term  is 


-b  ~  r  -us2 

2j  a_x  e  T  /  e 


ds 


2j  a 


-I3T?  /tT 
-1  6  T  VS 


11 


van  der  Waerden  shows  that  "the  remainder  of  the  integral  at  this  branch 

V- 

point  is  of  order  f 

■'  .m  - 3/2 

R  =  0(e  t  B  7  ) 

The  exponent  at  the  remaining  branch  points  are  larger  than  7J/T  and  the 
contribution  of  these  terms  may  be  neglected  for  B  large.  Thus,  we  have  the 
result 


-  e-w 

ku<b)  * v  T 


*/%  [(1-^)  cot2^  sinh 


ri+rn 


(22) 
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4,  THE. SOLUTION  OF  THE  SEMI-INFINITE  BIFILAR  HELIX 
4,1  The  Modified  Circuit  Equation  for  the  Semi-infinite  Helix 

The  semi-infinite  bifilar  helix  will  be  taken  as  two  helical  conducting 
wires  originating-  in  the  z"  =  0  plane,  of  a  cartesian  "coordinate  system  and  ' 
extending  in  the  positive  z  direction.  There  will  be  no  conducting  structure 
for  negative  values  of  the  z  coordinate,  and  the  structure  will  be  considered 
to  be  coincident  with  that  described  in  the  preceeding  section  for  positive 
values  of  the  z  fpoordinate .  The  lack  of  conducting  structure  in  the  nega¬ 
tive  half-space  forces  the  current  distribution,  I  (z),  to  be  identically 
zero  there.  Since  there  is  no  conducting  boundary  in  the  negative  half-space, 
the  electric  field  strength  tangent  to  the  null  line  in  this  space,  E  (z),  need 
no  longer  be  zero.  However,  E  (z)  must  be  identically  zero  for  positive  z. 

The  positive  subscript  is  used  to  indicate  a  function  of  z,  identically 
zero  for  all  negative  z.  The  negative  subscript  denotes  a  function  which 
vanishes  for  positive  z. 

Equation  (10),  in  this  case,  becomes 


V(z)  =  j  ?./k  f  G11(z>  z 0  sin  ^  ^  dz;/ 


dV  "  p 

sin  +  Z-  =  -E_(z)  -  jk?>  \  G11(ZJ  z 0  ^(z)  •  u^z')  1^  (z')  dzy 

.  -oo  + 

where  only  one  wire  is  considered.  The  potential  V(z)  is  eliminated  by 
multiplying  the  derivative  of  the  first  equation  by  sin  *1*  and  equating  to  the 
second  equation,  The  result,  after  a  simple  algebraic  manipulation  is 


oo 

-J 


-j  k/£  E_(z)  =  ’  |  Z^Cz-z')  *’I1  (zA)  dz' 

-00  ^  * 
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This  Is  the  modified  circuit  equation  for  the  semi-infinite  helix.  The 
circuit  equation  for  the  semi-inf inite" bifilar  helix  is  obtained  from  this  by 

‘  i!  _  "  .  •  •. 

including  the  contribution  of  the  second  wire  as  was  done  in  Chapter  3. 

"  "  ’i  ”  "  " 

Replacing  ,jthe  expression  on  the  left  with  V  (z)  giv^ps  the  result 


V  (z)  =  J  Z(z-z')  I  (z')  dz' 


for  the  circuit  equation  of  the  semi-infinite  bifilar  helix. 

4.2  The  Source  Problem 

The  Fourier  transform  of  the  convolution  integral  Equation  ($5)  is 

V  <B)  =  Z(S)  I  (3) 

1  t 

This  equation  contains  two  unknown  functions  of  3.  Since  the  support  of 

V- (z)  is  the  negative  z  axis,  its  transform  w 

co 

V-tB)  =  -i  J*  V_(z)  ejSzdz 
-eo 

is  regular  in  the  lower  half  of  the  complex  3-plane;  and  since  the  support 


of  I  (z)  is  the  positive  axis,  its  transform 
+ 


s  J  h 


<z)  eJBZdz 


is  regular  in  the  upper  half  of  the  complete  3- plane-.  It  will  be  assumed 


that  V  (z) j  and  I+(z)  are  of  decaying  exponential  order  at  infinity  so  that 

their  regions  of  regularity  include  the  real  B-axis .  The  procedure  for 

■* 

solving  Equation  (26)  used  in  this  study  depends  upon  factorizing,  Z(6)  in 


the  form 


z(S)-=  -TsT 


where  Z+ (3)  is  regular  in  the  upper  half-plane  and  Z  (B)  is  regular  in  the 

(12) 

lower  half-plane.  By  the  W.iener-Hopf  technique  f  Equation  (26)  can  be 
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written  >; 

V  (B)  Z  (6)  =•  I  (3)  Z  (B)  =  J ( 6)  W  (28)  ' 

-  t  .  t 

.  ■-*  ‘ 

If  the  factors  of  Z(G)  have  a  common  strip  of  regularity  including  the  real 

6-axis,  then  each  side  of  Equation  (28).  must  be  the  analytic  "continuation  of 
the  Other  in  its  region  of  regularity,  and  it  follows  that  each  side  is  an  entire 
function,  J(B).  If  this  function  has  polynomial  growth  at  infinity,  ..it  is'  a 
polynomial,  say  P^CG),  by  Liouville’s  theorem.  The  factorization  of  Z(B)  is 
not  unique  because,  Z^  (B)  and  Z  (6)  may  be  multiplied  by  the  same  polynomial, 
Q^(G),  without  changing  Equation  (27)  , 

To  obtain  a  unique  solution  of  Equation  (28)  it  is  necessary  to  determine 
the  asymptotic  behavior  of  the  factors  at  large  Values  of  B.  This  is  equi¬ 
valent  to  determining  the  behavior  of  J  (*)  near  the  origin.  In  this  problem, 
this  determines  the  source  conditions  on  the  current  distribution. 

The  backfire  bifilar  helical  antenna  will  be  fed  at  z  =  0  from  a  balanced 
transmission  line  on  its  axis.  In  the  mathematical  model  considered  here,  the 
transmission  line  and  the  wires  connecting  it  to  the  helix  will  be  neglected, 
and  the  current  dis l i ibu lion  I  (z)  will  start  from  some  finite  value  at  the  origin 
This  discontinuity  in  the  current  at  the  origin  implies  an  oscillating  point 
charge  distribution  at  the  end  of  the  wire.  An  infinitely  long  thin  straight 

wire  driven  by  an  oscillating  point  charge  at  its  end  has  been  discussed  by 
13 

Schelkuncff  .  The  point  charges  can  be  considered  as  the  source  of  the  currents 
flowing  in  the  helical  wires.  They  will  be  replaced  in  the  physical  antenna  by 
the  transmission  line  and  the  connecting  wires  . 

The  discontinuity  in  the  current  distribution  at  the  origin  implies  that 

r*  .  /aL-77. j,  !  i 

the  transform  1^,(6)  behaves  as  1/B  at  infinity.  The  current  di  stribut'ion 
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!«  ..  ,,  -  „  BQ  (P) 

...  x4-(p)  =  z  (P)  q  (p)  =  z rw> 

4  n  + 

Thus,  although  the  factorization  of  Z(Pj)  is  not  unique,  the  solution  for 

'V 

I  .(P),  based  upon  the  source  condition 

I  (z)->  I  (0) 

+  •+..- 

z - *  0  , 

is  unique. 

4,3  The  Wiener-Hopf  Factorization  ‘ 

12 

A  function  K(P)  can  be  factorized  by  theorem  C  of  Noble  in  the  form 

K(P>  =  K+(P)  K_(P) 

provided  that  K(P)  is  regular  and  non-zero  in  a  strip  containing  the  real 
P-axis  and  further  that  K(P)  tends  to  +1  as  the  magnitude  of  P  becomes  large 
on  the  axis.  The  required  factors  then  are  given  by  Cauchy's  integral  formula 
as 


oo 

K+(P)  =  exp  [-—■  j *  dxL  >  0 

-OO  .  o 

(31) 

oo  ' 

K_(P)  =  exp  [~j  jf  dx),  >m(P)  <  0 

-OO  * 

These  functions  are  bounded  and  non-zero  in  their  regions  of  regularity,  and 
they  tend  asymptotically  to  +1  along  the  real  axis,  This  theorem  will  be 
used  to  factorize  Z(p)  after  it  has  been  properly  conditioned. 

Some  typical  graphs  of  the  function 


47 T2 


(P)  =  t<7)2  "  (J)  ]  Z  B  (P)  -  ■  S  [  B  (P+T  )+B  (P-t)]  (32) 

t  ,  ,  n  n  n 


n  odd 


V*  \ 


n 
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for  real  values  of  3  are  plotted  in  Figures  8 -and  9.  In- Figure  8  the  frequency 

is  below  the  critical  frequency  Ck<kc),  and  in  Figure  9  the  frequency  is  above 

the  critical  frequency  (k>k  ).  In  the  first  instance  there  are  three,  roots  of 

c 

Z(B)  on  the  positive  13-axis,  and,  of  course,  three  corresponding  roots  on  the 

negative  8-axis  since  Z(B)  is  an  even  function  of  B.  The  analysis  given  here 

will  be  restricted  to  the  second  case,  where  the  frequency  is  above  the  critica 

value  k„,  and  where  there  is  only  one  real  root,  B  on  the  positive  B-axis. 

Y  o 

The  rdason  for  this  restriction  is  the  experimentally  observed  fact  that 

below  the  critical  frequency  a  current  traveling  wave  corresponding  to  the 

smallest  root  of  Z(8)  dominates  the  current  distribution.  This  result  is  shown 

in  Figure  11  of  Chapter  6.  The  dominance  of  the  contribution  of  this  root  is 

explained  by  considering  the  residue  of  the  inverse  of  Z(B)  at  these  roots. 

The  residues  are  inversely  proportional  to  the  derivatives  of  Z(B)  at  its  roots 

The  derivative  at  the  smallest  root  is  always  smaller  since  the  function  is 

logarithmic  in  the  vicinity  of  the  two  larger  roots,  and  therefore,  its  residue 

is  the  largest.  Above  the  critical  frequency  this  root  has  vanished  and  only 

the  relatively  weak  root  at  B  remains. 

o 

The  branch  points  of  Z(B)  at 

8  ■=  nT  4-  k  (n  =  »•  1,  +  3,  .  .  . ) 

all  lie  on  the  real  axis.  This  makes  it  impossible  to  .factorize  the 
function  Z(B)  into  factors  with  a  common  strip  of  regularity  including 
the  real  axis.  This  difficulty  can  be  surmounted  by  assuming  that  the 
medium  is  slightly  dissipative  "  ■ 


The  location  in  the  complex "B-plane  of  the  branch  points  and  zeros  of 
Z(B)  for  a  dissipative  medium  is  shown  in  Figure  Iff.  The  branch  cuts 

■■  ~  ..  -  \ 

..  .  '  ■'  ..  Y 

-  -  / 


Figure  8.  TUs  Bifilar  Helix  Determinant  a  1  Funct ion.  k  <1 


;ure  9.  The  Bifilar  Helix  Determinants!  Function  k  >  k 
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shown  are  chosen  to  satisfy  the  condition 

If  B  as  £  +j7]  the  equation  of  the  branch  line  is  „ 

(i-nT)TJ  =  -k a 

rv 

This  provides  a  strip  of  regularity  of  Z(B)  of  width  2a  including  the  real  axis. 
The  function  Z(8)  can  be  written  as  the  product  of  two  factors 

Z(6)  =  R(8)  K(S)  (33) 

such  that  K(B)  is  non-zero  in  the  strip  and  tends  to  +1  as  8  becomes  large. 

The  factor  K(B)  is  chosen  so  that  it  can  be  factored  by  Equation  (31) . 

It  is  particularly  chosen  so  that  its  function  values  are  real  over  as  much  of 
the  real  B-axis  as  possible.  This  allows  a  relatively  simple  estimate  for  the 
magnitude  of  K+(B)  to  be  obtained  as  shown  in  Equation  (41) ,  The  factor 
R(B)  remaining  will  be  factored  by  a  variety  of  devices.  These  factors  are 


R(B)  =  A 


b2-b2 

o 


(S2-62) 

o 


1/4 


(34) 


and 


K(B)  = 


(s2+k2)1/4 

(B-B‘ )  A 

O 


(35) 


r*  2  2 

B  is  the  root  of  Z(B)  and  the  factor  (8  -8  )  removes  the  zeros  from  K(B) 
o  o 

satisfying  one  of  the  requirements  of  Theorem  C.  The  exponential  factor 

12 

is  chosen  to  give  the  correct  asymptotic  form  and  to  be  easily  factorizable 


The  factor  (B2+k2)  is  chosen  to  give  asymptotic  behavior  like  /  B^2  J  . 

This  factor  is  chosen  because  it  is  real  on  the  real  B-axis.  The  term  k  in 
the  factor  is  arbitrary  as  any  other  constant  larger  than  a  would  make  this' fac¬ 
tor  regular  inside  the  strip  of  width  2a.  With  this  choice  of  factors  K(B)  is 
real  for  real  k  everywhere  on  the  real  axis  except  in  the  segments  |  B-nT|<  k. 


1/4 
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These  segments  will  be  called  the  "visible  regions"  of  the  real  (3-axis  since 
only  this  part  of  the  Fourier  spectrum  contributes  to  the  radiation  pattern  as 
is  shown  in  Equation  (4). 

The  functioh  K((3)  is  well  behaved  except  a  logarithmic  singularity  at 

I  (3  I  =  +  1  +  k 

This  singularity  is  integrable  in  Equation  (21),  and  K((3)  can  be  factorized 
by  Noble's  Theorem  C. 

4.4  Factorization  of  ■  tile 'Log- Integrable  Function 
If  K(|3)  is  represented  in  polar  form 

K((3)  =  M((3)  exp  [j(j>(p)] 

we  may  write 

..."  .  OO 

ve>  ■  afj  /  — ^  06) 

..  •  — OO 

The  integral  may  be  separated  into  its  real  and  imaginary  parts; 

°o  oo  M(x) 

/  ^  -■  -  -W>  *  /  d* 

-oo  M  -OO 

(37) 

+j(TT  In  M(p)  +  J  dx) 

-00 

The  imaginary  part  of  this  integral  gives  the  magnitude  of  K+((3), 

I  K+(f3)  I  »  V»(p)  exp  Zw  /  dxl  (38> 

-oo 

and  the  real  part  gives  the  phase.  The  computation  or  phase  involves  the 
evaluation  of  an  infinite  integral  while  the  magnitude  is  given  by  a. sum 
of  integrals  over  the  visible  range  only 
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oo  -■  >*  nT+k 

f  dX  =  e  r  4^5)  dx 

J  x-B  n=-co  y  x-B 

-oo  n>-k 


where  the  integral  over  the  range  containing  B  is  understood  in  the  sense  of 
the  cauchy  principal  value. 

To  show  that  the  contribution  to  the  phase  integral  of  all  terms  except 
for  n  1,  the  following  estimate  of  the  phase  is  obtained.  The  relative  size 
of  the  phase  values  for  n=l  and  n-2  can  be  obtained  from  Figures  8  and  9. 


The  function  Z(G)  may  be  written 


4ir  z(8)  =  E 

n  odd 


^n 


(B) 


where 


(B)  =  } 


(6  -k  ) I  (b  y  )K  (by  ) 

n  n  n  ‘n 


(39) 


If  B  is  in  the  n-th  visible  range  the  n-th  term  of  this  series  will  be  complex 


(B)  =  -  £<G2-k2)J  (b'y  )[N  (by  )+jJ  (by  )] 

n  z  nnnnnn 


where 


+7Tk2cot24' 


J  b'y  )[N  _ (by  )t  jj  .  (by  ) ] 

n  n  n-1  n  n-1  n 

I  J  <b/y  )  [N  (by  )+-jJ  _  (by  )  ] 

l  n+1  n  1  n+-l  'n  .  n+1  ’n 


r  ,2  T  3  ■,  1/2 

■y  =  [k  -(G-nT)  ] 


(40) 


For  large  values  of  n  the  first  term  dominates,  and  the  phase  of  is  given 


by 


<p  (B)  ~  tan 
n 


J  (by  ) 
-  n  !n 

N  (by  ) 
n  *n 


-  (<, 

( lB-nJ |  <  k) 


which  is  well  approximated  with  the  small  argument  form  of  the  Bessel,, functions 
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k  cot  41  <  <  nT 


The  result  is 


|»  <P)  ~  ~rr~TT, 

rn  n'.  (n-1)! 


[(?  - » 2-  (?) 2] 


J.  , 
cot  41 


which  tends  toward  zero  very  rapidly  when  n  increases.  The  remaining  terms 
of  the  series,  being  positive  real,  increase  the  real  part  of  the  function  in 
the  visible  range,  thus  making  the  phase  of  the  function  even  smaller  than  the 
phase  of  the.  complex  term. 

The  phase  of  K(P)  is  the  same  as  that  of  Z(P)  except  for  IP  I  <  P  ,  Here 
division  by  P  -Pq  makes  K(P)  positive  where  Z(P)  is  negative,  Only  in  the 
range  (P)  <  k  is  there  a  non-zero  phase  contributed  to  K(P)  by  the  exponential 
factor.  Therefore,  the  phase  integral  will  be  approximated  by 


/  jT 


n  =  -1,  0,  1 


is  in  the  range  (T-k,  T+k)  the  integral  in  that  range  is  written 


T+k 

C  Mx)  -  6(P) 

j  x-P 

T-k 


T+k 

^  /  ~{2 


>(x)  -  6(P) 


dx  +  ^(P)  in'  +  jIT6(P) 


Thus  the  magnitude  of  K+(P)  Is  given,  approximately  by 


T+k  +k  -T+k: 

f  ai  +  f  m  f  t 


-T+k : 

/  Kr1 


(41) 
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The  magnitude  of  (6)  is  given  by  two  integrals  over  a  finite  range 
while  the  computation  of  phase  involves  the  evaluation  of  an  infinite  integral. 
The  integrals  of  these  functions  are  not  known  in  closed/;|?orm,  and  they  must 
be  evaluated  by  numerical  tenchiques.  It  was  decided  to  evaluate  only  the 
magnitude  of  Z+(6)  for  this  reason,  and  also  because  fhis  is  sufficient  to 
determine  the  magnitude  of  the  radiation  pattern. 

A . 5  Factorization  of  the  Remainder  Function 

Siijtce  K+(8)  tends  toward  +1  for  large  8,  the  remainder  function 


2  „2  r 

nR(B)  =  A  -s - 2  1/d  exp  ( 

[B  +k ]  7 


B2-k2) 


must  be  factorized  to  give  the  desired  asymptotic  behavior 


R  (6)  ~  A6 

/l,  + 

l  ;• 

In  addition,  R^,  R  and  their  inverses  must  be  regular  in  the  indicated  half 
planes. 

The  zeros  (+13  )  of  R  can  be  assigned  to  the  appropriate  half  plane  by 


considering  the  medium  to  be  slightly  dissipative 


k  =  k  -j a 

■  ■■  o 

The  root  Bq  is  a  function  of  k  and  can  be  expanded  in  a  Taylor  series  about 

k  .  The  change  in  B  due  to  the  assumption  of  a  small  loss,  correct  to  the 
o  o 

first  order  in  a,  is  given  by 

9,3 

VJ??  *  Bo"ja  8T  . 

Since  B  increases  with  k,  as  seen  in  Figure  6,  the  root  at  +B  is  shifted 

..O  '  o 

to  the  lower  half  plane  by  a  small  loss,  and  thus,  belongs  to  R+.  Similarly, 

the  root  at  -B  is  a  root  of  R  . 

..  ■■  o  ....  - 

The  exponential  factor  .of  R  may  be  factorized  by  a  method  due  to 
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W  S  Ament  given  in  Noble 

By  this  method,  the  radical  in  the  experiment  is  written  in  the  form 

-VP2-k2  =  (k2-P2)  (P2-k2) 

where  the  first  factor  is  entire,  the  transformation  P  =  k  cos  0  is  used  to 

// 

remove  the  branch  points  of  the  second  factor 


g  =  (P2-k2) 


jk  sin  0 


As  shown  in  Figure  10,  the  branch  point  at  P=k  is  in  the  lower  half  plane  and, 
therefore,  must  be  assigned  to  g  .  This  branch  point  corresponds  to  the 
pole  at  0  =  0  in  the  0  plane.  Therefore, 

+  _  17-0  0 

g-g++g_  gin  Q  ^  sin  0 

is  the  required  separation  of  g.  When  these  functions  are  transformed  back 


to  the  P-plane, 


<P)  = 


77+ j  In 


p+Vp2- k2 


7T  V  P2-k2 


has  a  branch  point  only  at  P=k  and 


.  .  p+YP2-k2 

-J  In - - - 

g_(P)  =  - —4 - 

777  P-k2 

has  a  branch  point  in  the  finite  P-plane  only  at  P=-k.  The  factor  of  the 
exponential  function,  that  is  regular. in  the  upper  half  plane,  then  is 

-T|  Vp2-k2  p+Vp2-k2 

exp  17  +  j  In  < 

This  function, has  unit  magnitude  for  ail  negative  P  less  than  -k. 


The  remaining  factor  of  R,  can.be  written  P-jk 


The  pole  in  the  first  factor  is  in  the  upper  half  plane^  so  this  factor  is 
part  of  H  .  .The  contribution  of  the  second  factor  to  can  be  obtained  by 


Noble's  theorem  C  as  in  Section  4.4'.'  If  this  factor  is  called  G,  G  is 

+ 


given  by 


G^ (6)  =  exp 


(s  / 


and  its  magnitude  is  given  by  the  principal  value  of  this  integral. 

The  magnitude  of  the  remainder  function  then  is  the  product  of  Equation  (42) 


and  Equation  (43)  with  the  zero  at  13.  This  may  be  written  explicitly  as 


IB  (6)1  =  16-13  |  |  exp[-7)  Vf32-k2  (ir+j  In  8+V»2-k2)]  expf 
+  ° .  "  i r  "  k 


-  k  PJ  cot~I  i  dxJ 


-oo  x-6 


The  method  used  for  the  numerical  evaluation  of  Equations  (41)  and  (44)  is 


discussed  in  the  next  chapter. 


5,  EVALUATION  OF  THE  SOLUTION  BY  NUMERICAL  ANALYSIS 

“  ■  '  ’’  tl  •  If 

5.1  The  quadrature  Formula 

14 

A  Gaussian  quadrature  formula  is  used  in  this  study  to  evaluate  the 
integrals  in  Equations  (41)  and  (44).  It  is  used  because  it  offers  a  potential 
accuracy  comparable  to  other  available  quadrature  formulas  while  using  only" 
half  as  many  ordinates.  This  is  an  important  consideration  here  because  of  the 
complexity  of  the  integrand.  Ir  selecting  a  particular  Gaussian  quadrature 
formula  for  an  integral  of  the  form  ■ 

■  b  .  . " 

I  =  £  w(x)  f  (x)  dx 
a 

attention  must  be  given  to  the  weighting  function  P(x), 

For  the  weighting  function  w(x)  =  1  and  the  interval  (-1,  1),  the  Legendre 
Gauss  quadrature  is  used.  This  quadrature  formula  has  been  well  tabulated^ 
and  can  be  easily  instrumented.  The  error  function  for  this  quadrature 
formula  is 

E  _  2  (ml )  (2m) 

(2dh-1)[  (2m)'  ]3  f  (g) 

where  £  is  in  (-1,  1).  The  integrand  in  Equation  (41)  behaves  as 

1 

+(x)  “  In(l-x) 

as  x  approaches  -1,  The  first  derivative  of  this  function  and  all  higher 
derivatives  are  unbounded  near  the  end  points  of  the  interval.  Since  t,  can 
be  any  point  in  the  interval,  no  bound  on  the  error  of  the  quadrature  formula 
can  be  obtained.'  This  situation  could  be  corrected  by  using  the  weighting 
function  -  ■  y 
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in  a  small  region  near  the  endpoints  of  the  interval  (1-k,  1+k) .  « 

The  use  of  a  logarithmic  weighting  function  near  the  end  points  of  the 

■  ■  .  .  r  .  • 

interval  would  require  the  construction  of  a  new  set  of  orthogonal  polynomials 

On  the  interval  of  integration.  This  involves  the  moments  of  the  weighting 

a  f  " 

function  which  must|be  evaluated  in  terms  of  the  exponential-integral  functions. 


Since  this  difficulty  is  caused  by  the  behavior  ,0f  the  integrand  near  the  endpoints 
of  integration,  and  since  the  contribution  of  the  end  regions  to  the  integral  . 
is  obviously  small,  it  was  decided  to  use  a.  sixteen  point  Legendre, -Gauss  quadra¬ 
ture.  Although  no  bound  can  be  set  for  the  error  of  this  quadrature  formula 
It  is  reasonable  to  expect  that  the  error  will  be  small. 

5.2  The  Statement  of  the  Problem  in  the  Form  Used  by  the  Automatic  Digital  Computer 
In  preparing  a  problem  for  Automatic  digital  computation,  greatest  attention 
must  be  given  to  the  most  frequent  computations  in  the  program.  The  Legendre- 
Gauss  quadrature  of  the  integral  in  Equation  (41)  can  be  considered  as  controlling 
the  program  from  the  point  of  view  of  running  time.  This  is  caused  by  the 
comparatively  long  time  required  to  calculate  the  eighteen  Bessel  function 
values  required  for  each  ordinate  point.  Since  the  quadrature  points  do  not 
depend  on  the  particular  values  of  B  for  which  the  integral  is  evaluated,  it 
was  decided  to  compute  the  ordinate  values  in  a  separate  program  called  the 
"table  generator".  This  program  computes  the  function  in  Equation  (17)  subject 
to  the  approximation  in  Equation  (19)  at  each  of  the  sixteen  quadrature  points 
and  at  the  thirty-six  observation  points  given  by 


S  =  T+k  cos  6 


where  0  takes  on  increments  of  5°  from  5°  to  180°.  If  F(— )  =  Z(B)  is  used 

..  T2 

to  represent  the  approximate  form  of  the  bifilar  helix  determinantal  equation, then 


iff 

u 


k2  2a, 

„  „  .  —  cot.Y 

2  2  2 

><T>  -  2  =-[.<L  -  ;^>  8aC8)  -  T— -  ffl1.1«)tBwl(8)) 
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(45) 


B 

-x  — 
-e  T 


„2  ,  2  ,  -x  y3  -3x 

i/B  k  .  , ,  ,1-ye  .  „  ,  e  -x.  •, 

|(  -  — )  [in  (  ■■*—_-)  +2  (— - +  ye  )] 

t  t  i+ye 


—  cot2^  cosh  X  [ln(l-y2e_2x)  +  — % -  +  yV2*] 

t2  a 


Vsin..^  sin  +  ' 


with  i  =  +  1,  +  3 


is  the  expression  used  in  the  machine  computations,  where 


_  1 _ 1 

X  “  sin  41  sin  4|/ 

_  cos  ^  (1+sin  41) 
y  ~  cos  4*  (1+sin  41') 

k  P 

The  input  data  for  the  table  generator  is  the  normalized  frequency  (;p  =  , 

the  cotangent  of  the  pitch  angel  4'  and  the  relative  wire  size  6.  The  constants 
x  and  y  in  Equation  (45)  are  then  computed  as  are  those  functions  of  these 
constants  not  depending  on  13,  These  are  stored  for  use  in  the  later  calculations. 
The  eight  roots  of  the  sixteenth  Legendre  polynomial  in  (0,1)  are  then  used  to 
compute  the  values  of  F(B/T)  in  the  range  (0,  l~k/T) •  The  results  of  this 
calculation  together  with  the  sixteen  .quadrature  points  and  thirty-six 
observation  points  in  the  range  (l-k/T,  l+k/T)  are  emitted  on  punch  cards  for 
use  by  the  subsequent  integrating  program.  This  program  also  finds  the  zero, 

B  3  ' 

o,  of  F(6/T)  that  must  be  used  in  the  subsequent  calculation.  .. 

T 

The  zero  seeking  routine  used, in  this  program  starts  with  an  initial  value 
of  the  argument  increment  size,  and  direction  of  search.  The  increment  size  is 
doubled  at  each  step  until  the  function  values  change  sign.  When  this  occurs  a 

second  approximation  to  the  root  is  obtained.  The  value  of  the  function  at  this 

. 

approximate  root  is  calculated,  and  ff  it  is  larger  than  10  ,  the  search  is  again 


Initiated  toward  the  trug  root  location  with  an  increment  one  tenth  as  large 

.  ’  "  J.V 

as  the  last  increment  used.  This  process  continues  until  either  the  function 

—  g  O 

value  at  the  approximate  root  is  sufficiently  small  (less  than  10  ).,  or  until 

the  increment  "becomes  too  small  to  modify  the  argument  of  the  function  (less 
than  lone  in  the  eighth  significant  figure), 

'  M 

The  second  and  final  program  in  the  set  computes  I+(T+k  cos  0)  and 

I+(~T+k  cos  0)  using  the  tabulated  data  provided  by  the  table  generator. 

"  ,  O' 

For  each  value  of  0  in  5  increments,  it.pomputes  the  following  factors. 


P  (0) 
o 


J5._  ?o 


P  (0)  = 
o 


.  p  <  k 


P  >  k 


P2(0)  =  e 


J 


p  cot  1  (|) 


P3(0)  = 


[<e/T>W)2]1/4  e  W(?)2-(f)2 

(pa)2-(|^)2  »t  t 

i+— 


P.(0)  =  e 
4 


l+f  <*> 

£/I  I  2 

IT  J,  k  X  -(S/T) 


where 


0  ,  k 

=  1+—  cos  0 
T  T 


The  constant  0  in  the  first  factor  is  the  one  positive  real  root  of  Z  (0). 

o  + 

and  is  supplied  by  the  table  generator.  The  variable  0  is  generated  in  the 

o 

Integrating  program  at  5  increments  in  the  polar  angle  0  corresponding  to  di¬ 
rection  of  observation  of  the  antenna.  The  first  factor  can  be  interpreted 
as  the  contribution  of  the  higher  waveguide  mode,  whose  propagation  constant 
is  Po,  to  the  radiation  pattern.  The  next  two  factors,  P^G)  and  P2(0),  are 


included  to  insure  the  proper  asymptotic  behavior  of  Z  (B)  and  consequently  of  .. 
I  (6) -  The  values  are  generated  in  the  integrating  program  independently  of  the 

+  B  '  ' 

function  values  of  F  (^)  supplied  by  the  table  generator.  These  factors  have 
relatively  little  effect  on  the  value  of  I+(B)  in  the  visible  range. 

The  last  two  factors,  P  (0)  and  (01.  exercise  primary  control  over  the 


shape  of  the  radiation  pattern.  The  fac^r  P  (0)  depends  primarily  upon  the 
magnitude  of  F(^)  at  the  observation  points.  P^C©)  depends  upon  the  phase  at 
the  observation  points  as  well  as  the  phase  at  the  quadrature  points.  The 


IV 


integrating  routine  obtains  the  real  and -imaginary  parts  of  F(— )  from  the  table 
supplied  by  the  table  generator  and  converts  them  to  polar  form.  The  magnitude 
is  then  used  to  compute  P^  while  the  phase  is  integrated  to  obtain  the  value  of 

P4* 

The  product1, of  these  factors  is  taken  for  the  magnitude  of  Z+(B)  .  The 
inverse  of  the  product  is  then  interpreted  as  the  magnitude  of  the  Fourier 
spectrum  of  the  current  distribution  The  vector  potential  for  the  backfire 
bifilar  helix  results  when  the  current  spectrum  is  multiplied  by  the  function 

J  cot  41  sin  8) 

as  indicated  in  Equation  (5).  The  radiation  pattern  can  be  deduced  from  this 
as  indicated  in  Chapter  2.  Radiation  patterns  were  computed  for  a  range  of 
pitch  angles,  relative  frequency  k/T,  and  relative  wire  sizes  The  results 

of  these  computations  are  discussed  in  the  next  section  where  they  are  compared 

-  .  X 

with  experimental  results. 
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-  r»  6.  experimental:  and  theoretical  results 

6,1  The  Experimental  Study  , 

The  study  of  the  bifilar  backfire  helical  antenna  was  started  as  part  „ 

ik 

of  a  larger  investigation  of  the  properties  of  periodic  radiating  structures, 

This  investigation  in  turn  was  started  in-  an  effort  to  learn  more  about  the  -• 

10_18  t 

behavior  of  log-periodic  antennas  .  In  this  connection  it  is  suggested 

that  the  properties  of  a  log-periodic  structure,  as  a  function  of  distance 

from  the  apex,  are  related  to  those  of  a  periodic  structure  whose  period  is 

given  by  the  local  period  of  the  log-periodic  structure.  In  this  sense  the 

18 

bifilar  helix  is  an  analog  of  the  two  arm  equiangular  spiral  antenna 

The  portion  of  a  log-periodic  antenna,  nearest  the  feed  point  at  or  near 
the  apex  of  the  structure,  acts  as  a  transmission  line  carrying  the  energy  to 
the  larger  portion.  The  energy  is  carried  to  the  so-called  "active  region" 
of  the  structure  whose  position  and  size  varies  ..linearly  with  frequency.  This 
region  is  thought  to  be  primarily  responsible  for  the  radiation  from  the  antenna. 
Beyond  this  region,  the  current  decays  rapidly.  These  characteristics  of  the  " 
log-periodic  antenna  are  observed  for  the  bifilar  helix  when  the  variation  of 
distance  from  the  apex  is  replaced  with  a  variation  of  frequency. 

The  operation  of  the  helix  as  a  waveguide  is  discussed  in  Chapter  3. 

There  it  is  seen  that  the  propagation  constant  is  given  approximately  by 

S  =  k/sin  + 

until  the  edge  of  the  visible  range  is  approached.  That  is  until 

„  » 

— =  ~-k 
sin  v 

"  n 

This  equation  defines  the  cutoff  frequency  or  critical  frequency  of  the 
principal  waveguide  mode  of  an  infinitely  thin  helical  conductor.  This  critical 


'  "ii 


frequency,  normalized  with  respect  to  is  called^C  to  avoid  contusion1  with 
the  critical  frequency  for  a  finite  size  helical,  conductor.  It  is  given., 
explicitly  by 


x  = 


sin  4* 


(46) 


1+sin  Y 

The  critical  frequency  marks  the  boundary  between  the  frequencies  for  which 
the  helical  structure  is  primarily  a  waveguide  and  the  frequencies  for  which 
it  is  a  useful  antenna. 

Several  models  of ' the  bif liar  helical  antenna  were  constructed  for  a 
range  of  yC.  between  0.05  and  0.4.  The  physical  dimensions  of  these  models 
are  given  in  Table  I.  These  are  based  upon  a  frequency  of  1.5  Gc ,  The 
pitch  of  the  helix  is  determined  from  the  wavelength  at  the  base  frequency  by  " 

P=JA  (47) 


The  pitch  angle  is  given  by 

4*  =  sin-1(^3^-) 

and  the  mean  radius  of  the  helix  is  given  by 


(48) 


(49) 


The  waveguide  operation  of  these  models  was  studied  by  sampling  the  fie'lds 
near  the  antenna  with  a  current  loop  moving  pa-rallel  with  the  axis  of • the  helix. 

A  typical  amplitude  response  of  this  loop  as  a  function  of  distance  from  the 
feed  point  is  shown  in  Figure  11.  The  antenna  was  fed  by  a  balanced  two  wire 
line.  Below  the  critical  frequency,  the  propagation  constant  of  the  waveguide 

...  ’  V.  |(  •  "  "  •• 

mode  can  be  easily  determined  by  the  standing  wave  on  the  structure.  The  propa¬ 
gation  constants  measured  in  this  way  are  sKown  in  Figure  12.  It  is  noted  that 
the  cutoff  frequency  for  each  of  these  models  is  somewhat  below due  to  the 
finite  wire  thickness,  The  cutoff  frequency  for  the  models  tested  in  this  invest! 
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Figure  12.  Measured  Propagation  Constants  oi  the  Principal 
Waveguide  Mode  on  the  Bifilar  Helix 
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Table  I 


The  Physical  Parameters  of  the  Experimental  Models  of  the  "Blf liar  Helical  Antenna 


Code  Name 

X 

Pitch  p(cm) 

Radius  b(cm) 

Relative* 
Wire  Size  6 

Pitch  Angle 
41  (degrees) 

BBH-1 

0.1 

2 

2.663 

.0227 

O 

.  6.38 

BBH-2 

0.15 

,3 

2.847 

„  .0242  1  ' 

10.18 

BBH-3 

0.2 

4 

2.465 

.0262  •' 

14.49 

BBH-4 

0.25 

5 

2.251 

.0287 

19.47 

BBH-5 " 

0.30 

6 

2.013 

.0323 

25.37 

BBH-6  „ 

0.35 

7  . 

1 . 743 

.0371 

32.60 

BBH-7+ 

0.40 

16 

2.85 

.0226 

41.82 

♦No .  16  AWG 

Tinned 

copperwire  was 

used  in  the  construction  of  the 

models . 

*  -(The  base  frequency  for  this  model  was  changed  from  1.5  Gc  to  .75  Gc  for 
mechanical  reasons. 


gation  is  plotted  as  a  function  of  the  pitch  angle  41  in  Figure  13,  where  it  is 

compared  with^,  the  normalized  cutoff  frequency  for  an  infinitely  thin  wire 

■  °  "  1/.:  -  " 

model.  n  ..  . 

Ati  unusual  feature  of  the  standing  wave  pattern  of  the  waveguide  mode  on 
the  bifilar  helix  was  noted  during  the  course  of  the  measurements.  The  apparent 
guide  wavelength  increases  with  frequency  from  a  minimum  at  zero  frequency  equal 
to  the  pitch  of  the  helix  to  a  maximum  at  the  cutoff  frequency,  the  maximum 
being  somewhat  less  than  the  free  space  wave  length.  This  indicates  that  the 
principal  waveguide  mode  on  the  bifilar  helix  is  a  backward  wave.  The  phase 
constant  decreases  with  increasing  frequency,  and  therefore,  the  phase  velocity 
and  the  group  velocity  must  have  opposite  signs.  The  group  velocity  must  be 
directed  away  from  the  feed  point  (i.e.  in  the  positive  direction),  therefore  the 
phase  velocity  is  in  the  negative  direction  as  indicated  in  Figure  12". 

The  measured  results  in  Figure  12  do  not  seem  to  agree  with  the  calculated 
results  show, n 'in  Figure  6.  However,  Figure  6  gives  the  phase  progression  of  the 
current  along  the  helical  conductor,  while  Figure  12  gives  the  phase  progression 
of  the  near  fields  of  the  bifilar  helix.  The  results  of  Figure  6  agree  with  the 

results  of  Marsh^  for  phase  measurements  along  the  wire.  When  measurements 

'  11 

along  the  helix  parallel  »fto  the  axis  are  made,  an  additional  phase  ^shift  of  2ir 
radians  per  turn  is  observed.  This  is  the  expected  result  since  the  test 
probe  is  now  moving  from  one  conductor  to  the  other  as  it  travels  along  the 
helix.  At  zero  frequency  therefore,  the  phase  changes  by  T T  radians  as  each 
conductor  is  crossed.  From  another  point  of  view,  the  additional  21T  radians 
of  phase  shift  results  from  the  rotation  of  a  tangent  to  the  helical  conductor 
as  it  moves  along  a  turn  of  the  helix. 

It  is  seen  in  Figure'll,  that  above  the  cutoff  frequency,  the  standing  wave 


has  vanished  and  is  replaced  'by  a  decaying  amplitude  distribution.  As  frequency 
is  increased  further  the  rate  of  decay  increases.  This  is  shown  in  greater  detail 
in  Figure  14.  At  these  higher  frequencies,  a  standing  wave  distribution  again 
appears  with  a  maximum  amplitude  about  twenty  decibels  below  the  input  level.  This 
is  the  higher  order  helical  waveguide  mode.  The  excitation  of  this  mode  increases 
with  frequency  from  a  negligible  level  near  the  cutoff  frequency  of  the  principal 
mode.  The  development  of  this  mode  places  an  upper  limit  on  the  frequencies  for 
which  the  bifilar  helix  behaves  as  a  backfire  antenna. 

Phase  measurements  made  in  the  near  field  above  the  cutoff  frequency  are 
shown  in  Figure  15.  In  this  figure  the  phase  is  normalized  to  the  free' space 
propagation  constant  and  is  plotted  in  centimeters.  The  first  curve  shows  that 
the.  phase  is  leading  as  the  probe  moves  away  from  the  feedpoint.  The  first 
curve,  slightly  above  the  cutoff  frequency,  k^, "shows  leading  phase  over  most  of 
the  length  of  the  structure.  As  frequency  is  increased,  the  direction  of  phase 
progression  near  the  end  of  the  structure  changes  and  is  away  from  the  feedpoint. 
The  point  at  which  the  phase  progression  changes  from  leading  to  lagging  moves 
near  to  the  feed  point  as  frequency  is  increased.  This  corresponds  to  the 
point  at  which  the  feed  region  currents  have  decayed  to  the  level  of  the  higher 
order  helical  waveguide  mode. 

The  near  field  phase  measurements  indicate  that,  although  they  do  not 
correspond  to  a  proper  mode,  the  feed  region  currents  have  a  backward  wave 
character  above  the  critical  frequency,  k^i  This  leading  phase  characteristic 
persists  as  "long  as  the  feed  region  currents  are  dominant  on  the"  structure . 

On  the  remainder  of  structure  the  higher  order  waveguide  mode,'  with  its 
lagging  phase  or  forward  wave  characteristic,  dominates  the  current  distri¬ 
bution.  The  fact  that  in  the  feed  region  the  direction  of  ohase  progression 


Figure  14.  Current  Amplitude  Distribution  above  the  Critical 


DISTANCE  fROM  FEED  POINT  IN  cm 


is  toward  th&  feedpomt^ 
bifilar  helical  antenna, 
the  pha.=e  progression  of 


is  consistent  with  the  "backfire"  radiation  of  the 
since  an  antenna  tends  to  radiate  in  the  direction  of 
current. 


The  measured  radiation  patterns  for  the  backfire  bifilar  helical  antenna 
are  presented  in  Figure  16.  The  patterns  are  plotted  so  that  their  centers  » 

give  the  helix  radius  and  pitch  in  wavelengths  at  the  frequency  for  which  the 

'  ■  {  '  .  • 

pattern  was  measured.  Two  curves  are  also  provided  for  reference'itq  the 

\  \  .. 

Brillouin,  diagram  for 'the  bifilar  helix.  The  curve  .  „ 

«  bA  y.iz'i.P'}. 

f  -  2  77 

gives  the  frequency  at  which  the  propagation  constant  asymptote 

S  ~  k/s  :in  41 


intersects  the  edge  of  the  visible  range  given  by 

S  =  T-k 


This  corresponds  to  the  cutoff  frequency  of  the  principal  mode  on  a  helical 
waveguide  with  infinitely  thin  conductors  The  curve 


b  V^l_t2_p/^ 

*  •"  2Tt 

gives  the  frequency  at  which  the  propagation  constant  asymptote  intersects 
the  edge  of  the  visible  range  given  by 

^  S  -  ti-k 


This  curve  is  given  by  Kraus  as  the  upper  frequency  limit  for  beam  mode 
operation  of  the  helical  antenna  The  range  of  parameters  for  this  antenna 
given  by  Kraus  .is  indicated  by  the  broken  line 

The  critical  frequency  can  be  determined  by  studying  the  radiation 
patterns  shown  in  Figure  16,  As  t.he  critical  frequency  is  approached  from 
above,  the  patterns  show  increasing  directivity  characterized,  in  general,  by 


a  single)  lobe  with  little  or  no  "etjd&Lre"  radiation  'or  back  lobe.  Below  the 


critical  frequency,  the  patterns  have  many  lobes  and  the  pattern  shape  changes 
rapidly  with  frequency.  "This  indicates  the  standing  wave  current  distribution 
on  the  antenna.  As  the  frequency  is  increased  the  patterns  tend  to  broaden  and, 
for  larger  pitch  angles!,  tend  to  form  a  split  beam.  This  is  due  td'  two  effects: 
first,  the  phase  progression  in  the  feed  region  tends  to  become  faster  than  light, 
and  second,  the  rate  of  current  decay  in  the  feed  region  increases,  decreasing 
the  effective  aperture..  The  splitting  of  the  main  beam  is  not  observed  for 
smaller  pitch  angles  because  the  effective  aperture  is  too  small  to  form1 distinct 
beams.  For  the  larger  pitch  angles,  the  rate  of  current  decay  is  smaller,  and 
the  larger  effective  aperture  can  produce  a  well-defined  beam,  This  result  is 
consistent  with  Equation  (20)  where  it  was  shown  that  the  coupling  of  the 
antenna  current  to  the  radiated  field  is  controlled  by  &/sin  + .  Since  the 
relative  wire  size  6  does  not  change  appreciably  over  the  range  of  models  Rested, 
the  decreasing  pitch  angle  serves  to  increase  the  coupling  and  hence  the  rate  of 
current  decay.  f,  »> 

6.2 _ Comparison  ofComputedPatterns  with  Measured  Patterns 

n 

The  effect  of  wire  size  is  not  shown  in  the  experimental  results  of  this 

study,  but  it  is  illustrated  in  the  computed  patterns  shown  in  Figures  17-23. 

II 

The  patterns  for  the  largest  wire  size  =  0.2)  show  the  least  variation  with 

1  .  ;  ,  '  •“  '  '  .. 

frequency.  The  patterns  for  the  smallest  wire  size  (°  =  0.002)  have  the; greatest  // 

"ii  ... 

directivity  near  the  critical  frequency,  but  tend  to  broaden  faster  with  frequency. 

It  chn  also  be  seen  from  the  computed  results  that  a  change  in  wire  size  is  more^-  v 

\\  \ 

effective  in  changing  the  shape  of  the  radiation  pattern  for  larger  pitch  angles. 


This  again  indicates  the  dependence  of  the  performance  of  the  backfire  bifilar 

c  "  |  ,  ij  i  i 

helical  antenna  on  °/sin  4* . 


Computed  Radiation  Patterns  for  BBH-1 
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Figure  18.  Computed  Radiattbn  Patterns  for  BBH-2 


K/t  *  2400  f=  1800  Me 

Figure  19.  Computed  Radiation  Patterns  for  BBH-3 


k/r  =  .3200  f  =  1600  Me 


k/r  =.3600  f  =  1800  Me 


k/T  =  .3400  f  =  1700  Me 

BBH  -  5 

8  =.2 - - 

8  =.02  - 

8  » .002  - - 


Figure  21.  Computed  Radiation  Patterns /for  BBH-5 


Figure  22.  Computed  Radiation  Patterns  for  BBH-6 


L 


frequency 


On'fe  pattern  was  computed  for  each  of  the  antenna  models  for  the  normalized 


=  % 


This  corresponds  to  the  model  frequency' of  1.5  Gc.  One  pattern  was  computed 
for  a  frequency  one  hundred  megacycles  below  this  frequency ,• and  three  were 
comput,,ed  in  increments  of  one  hundred  megacycles  above  this  frequency.  "  It 
was  considered  desirable  for  uniformity  to  use  the  same  wire  sizes  for  the  entire 


range  0  antenna  parameters.  The  value  (6  =  0.02)  was  chosen  as  a  convenient 

■1>‘  <  \  •  .  \\ 

number  near  tlie  values  obtained  for  the  experimental  models.  Wire  sizes  of 

one  tenth  and  ten  times  this  value  were  used  to  show  the  effect  of  wire  size 
\v  “  . 

on. the  performance  of  the  antenna.  ' 

The  computed  patterns  for  the  wire  s^ze  5  =  (0.02)  are  compared  with 

measured  patterns  in  Figures  24-30..  The  agreement  between  the  measured  pattern 

and  the  computed  patterns  is  good  in  the  main  beam.  There  is  a  considerable 

difference  in  the  shaped  of  the  computed  and  measured  patterns  in  the  direction 

of  the  structure  (back1 lobe) .  The  radiation  in  this  direction  is  controlled 

by  the  length  of  the  helical  antenna  in  the  manner  described  in  the  next 

chapter.  In  this  chapter  we  will  limit  our  discussion  to  a  comparison 

.  .  M  .  " 

betwe^eji  the  main  beams  of  the  computed  and  measured  patterns. 

In  the  patterns  for  BBH-1,  Figure  24,  the  agreement  between  the  computed 
and  measured  patterns  is  quite  good.  The  relative  wire  size  of  the  experimental 
model  6  =  0.0227.  This  does  not  differ  much  from  the  value  (6  =  0. 02)' used, in 
the  computation.  In  addition  it  can  be  seen  in  Figure  18  that  the  shape  of 
the  radiation  is  little  dependent  on  the  relative  wire  size.  For  BBH-2  through 
BBH-6  the  relative  wire  size  increses,  and  the  agreement  between  the  measured 
pattern  and  the  patterns  computed  for  5  -  o.02  becomes  poorer.  The  agreement  is 
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Figure  25.  Comparison  between  Measured  apd  Computed 
Patterns  for  BBH-2 


if 


Figure  30,  Comparison  between  Measured  and' .Computed 
Patterns  for  BBH-7 


best  for  the  lower  frequencies  where  the  relative  wire  size  has  less  control 
over  the  pattern  shape.,  Tfyjp  measured  and  computed  patterns  show  the  greatest 
difference  for  BBH-6  where  the  relative  wire  size  (&  0.037)  has  its  largest 

'  ,5>  •  . 

value.  When  BBH-7  was  constructed  it;  was  found  that  the  polystyrene  foam  core 
used  to  support  the  helical  wires  had  insufficient,  strength  when  sized  for  the 
design  frequency  of  1,5  Gc  Tile  design  frequency  was  changed  to  0.75  Gc, 

This  doubled  the  radius  of  the  helix  and  halved  the  relative  wire  size  to 
0  =  0.0226.  This  is  evident,  in  Figure  30  from  the  improved  agreement  between 
the  measured  and  computed  patterns  '  ■ 

Computed  values  of  the  axial  ratio  on  axis  in  the  main  beam  are  shown  in 
Figufe  "31.  These  are  obtained  from  Equation  (6)  (l 


AR 


1+r 

T-r 


// 


where  r  is  the  ratio  of  the  righr -handed  and  left-handed  circularly  polarized 
components  of  the  radiation  pattern  given  by 


1 

\ 


jo- y  . 


It-l-k.) 


u  -  «  • 

The  axial  ratio  is  best  (nearest  unity)  at  the  -^lowest  operating  frequency. 

These  results  are  comparable  with  those  obtained  with  the  helical  beam  antenna. 

The  comparison  of  the  measured  patterns  and  the  computed  patterns,  allow¬ 
ing  for  the  differences  in  relative,  wire  size,  shows  that  the  mam  beam  radiation 
.  -  „ ,  '  ..  11  , 
pattern  is  predicted  by  the  calculations  outlined  in  Chapter  5  based  on  the 

theory  of  Chapter  4.  This  theory  is  bNased  upon  an  infinitely  long  antenna.  The 
effect- of  the  length  of  the  structure  accounts  for  the  difference  between  measured 
and  computed  results  as  dtescribed  .  i n  the  next,  section. 


6.3  Endft.re  Radiation  from  the  Backf  ire  ..Bif  liar  Helical  Antenna 

The  endfire  radiation  from  the  backfire  bifiiap  helical  antenna  is 


Helical  Antenna 


controlled  by  the  root  (3  of  the  determtnantal  equation.  Values  of  this  root 

o  " 

for  the  parameters  used  in  the  pattern  calculations  are  shown  in  Figure  32. 

3  '  3 

In  this  figure  the  normalized  values  of  the  root.  |j=  ^2-1  have  been 

shifted  to  the  right  one  unit  to  account  for  the  difference  between  the 

apparent  phase  constant  as  one  travels  along  the  wire  and  as  one  travels  along 

'  \\.  ... 

the  helical  axis.  Thus  zero  corresponds  to  a  circularly  polarized  traveling 
wave  of  current  traveling  along  the  positive  z-axis  from  the  origin. 


I  e“(jBo  +  Q)Z 
b  o 


z  >  0 


The  pattern  produced  by  this  distribution  is 


P(0)  =  I 


- [B  -k  cos  0)+a]z  . 


13^-k  cos  6-j a 


For  vanishing  a,  the  magnitude  of  P(0)  at  9=0  (the  endfire  direction)  is 


(o)  B  *-k 
o 


\  .. 

In  this  equation  I  is  proportional  to  the  residue  of  I  (B)  at\\B  =  B  .  This  • 
o  +  \\  o 

\\  /x/  f/  ,, 

residue  is  inversely  proportional  to  the  slope  of  Z(B)  at  B  =  Bq.  Figure  10 

~ 

shows  that  Z(B)  has  a  logarithmic  singularity  at  B  =\+k.  When  B  is  close  to 

>y  ,  ^  V  ' 

T+k,  the  case  for  k  near  k  ,  the  slope  of  Z(8)  at  B  is  very  large,  and  I  is 

C  „  °  . ?  - 

small.  As  k  is  increased  the  slope  of  Z(B)  at  B  decreases  ariSif.'A  increases, 

o  .  //%  o 

■■  .  •  f  _x  .  V  ■ 

At  the  same  time,  (BQ-k)  is  very  large  when  k  is ' near"  k  and  decreases  rapidly 
as  k  increases.  The  computed  patterns  for  BBH-1  and  BBH-2  show  no  radiation  due 
.to  because  (B^-k)  is  too  large.  At  the  othlr  end  of  the  range  of  parameters,. 


Figure  32.  Calculated  Values  of  the  Propagation  Constant 
of  the  Higher  Order  Waveguide  Mode  on  the 
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BBH-7  shows  little  radiation  due  to  Bq  because  X^^is  small.  The  effect  of 
relative  wire  size  on  this  phenomenon  is  shown  ;/in  the  patterns  for  BBH-3,  ^ 
Figure  19.  in.  this  figure  the  curve  for  the  largest  wire  sj/fce  0  =0.2)  does 
not  have, .radiation  due  to  B  while  the  smallest  wire  0 . 002).  does.  Endfire 

'•  Q.  i;  .i  •>  .  <  .  '  ' 

tf-  ■  -  '  4  \\ 

radiation  due  to  B  is  also  evident  In  the  patterns  for  the, /bt^er  models. 

If  the  bifilar  helical  antenna  is  long  enough  to  support  the  feed 
region  current,  the  shape  of  the  main  beam  is’  not  affec-tbd  by  truncation  of 

(I  '  " 

the  structure.  The  radiation  due  to  Bq  is,  however,  controlled  by  the 
length,  L,  of  fhe  antenna.  Considering  only  the  incident  wave,  Equation  (50) 
must  be  modified  for  the  finite  antenna  as  follows 


Ve>  -  J  e 

o 


V  -[j(8/ -k  cos  6)+a]z. 


-[ j(B  -k  cos  6)+a]L 

'■  ‘t  T  e  ° _ ~l 

o  -[ j(8^-k  cos  0)+c] 

For  vanishing,  a,  the  magnitude  of  P  (0)  at  0  =  0.  i& 

« \  L  , 

Bin 

PL«»  -  ^  '.V!  -y.i,  ° 
o 

Thus  when(BQ-k)  L/2  is  small  the  radiation  ^n  the  endfire  direction  is  reduced 

by  the  factor  (B  -k)L„  \\ 
o  " 

PT(0)  =  (B  -k)L  P(0)  (52) 

Lt  O 


The  magnitude  of  P  (0)  at  other  angles  is  given  by 

Li 


P(0)  =■•  I  L  -2^- 

L  o  x 


where 


x  =  (B^-k  cos  0)  L/2 
o 
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sv 
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This  term  and  a  similar  one  due  to  the  reflection  of  the  higher  order  waveguide 

mode  from  the  end  of  the  antenna  account  for  the  back  lobes  in  the  measured 


patterns  and  "scallops"  observed  on  the  main  beam  at  the  higher  f requencied'. 

.4  ..  }■  "I  "  -  ■  •  „ 


No, .attempt  is  made  here  to  obtain  an  exact  theoretical  preduction  of  the  effect 
of  the  length  of  the  antenna  since  this  requires  a  knowledge  of  the  far  field 
phase  pattern  of  the  feed  region  currents. 

The  effect  of  helix  length  on  the  performance  of  the  bifilar  helical  antenna 
was.  studied  experimentally.  Figure  33  djpws  some  ^ypical  patterns  resulting 
from  changing  the  length  of  the  backfire  bifilar  helical  antenna.  The  size  of 
the  back  lobe  decreases  with  length  with  little  change  in  the  main  beam  until 
the  antenna  becomes  too  short  to  support  the  feed  region  currents.  Shortening 
the  antenna  beyond  this  point  will  generally  broaden  the  main  beam  pattern. 

"  -H 

Front-to-back  ratios  obtained  for  different  lengths  are  shown  in  Figure  34. 

Figures  33  and  34  are  for  BBH-2  at  a  frequency  of  1.4  Gc.  A  similar  set 

of-  measurements  was  made,  at  the  other  end-  of  the  parameter  range,  for  BBH-7. 

For  this  antenna  there  was  no  measurable  endfire  radiation  for  any  of  the 

lengths  used  in  the  measurements.  This  result  is  explainable  when  the  calculated 

value  of  (3  is  observed.  For  BBH-7  at  k  =  o.4T 
o 

S  - (1+k)  <  lxlO-7 
o 

it  would  take  a  very  long  antenna  indeed  to  show  endfire  radiation.  In 
general,  a  greater  length  of  antenna  is  required  at  the  lower  frequencies  than 
at  the  higher  frequencies,  because  of  the  slower  rate  of  current  decay, 
observed  in  the  near  field  amplitude  measurements,  near  the  critical  frequency. 

At  frequencies  well  above  the  critical  frequency,  there  is  an  optimum 
length  for  the  bifilar  helical  antenna  for  each  frequency.  This  optimum 
length  is  just  that  required  to  support  the  feed  region  currents.  A  longer 
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Figure  33.  Radiation  Patterns  Showing  the  Effect  of 
Changing  the  Length  of  the  Antenna 


length  will  decrease  the  front-to-back  ratio  by  ^increasing  the  radiation 

due  to  8  .  At  frequencies  near  the  critical  frequency  the  feed  region  current 

o  ..  .. 

penetrates  the  structure  to  a  much  greater  depirh,  The  lowest  frequency  that  can 
be  used  in  the  backfire  mode  on  a  particular  model  depends  upon  the  length  of 
the  antenna.  The  required  length  changes  very  rapidly  as  the  critical  frequency 
is  approached  so  that  the  cutoff  frequency,  determined  from  the  pattern  measure¬ 
ments  on  an  antenna  of  a  reasonable  length  £10-15  turns) ,  is  in  error  by  less 
than  one  percent.  An  infinitely  long  antenna  is  required  to  support  the  feed 
region  current  at  the  cutoff  frequency. 

6.4  Other  Experimental  Results 

Pour  investigations ,  supplemental  to  the  main  body  of  this  study,  are' 

reported  in  this  section.  These  investigations, ,a:re  by  no  means  exhaustive 

1  •  '*  // 

and  serve  to  suggest  further  work.  The  first  of  these  is  a  study  of  the  effect 
of  placing  a  conducting  core  on  the  axis  of  the  helical  antenna.  The  second 
investigation  is  a  study  of  high-grain  backfire  bifilar  helical  antennas. 

In  the  third  experiment.,  a  long  slowly  tapered  helix  is  studied,  and  the  last 
study  shows  that  a  monofilar  can  be  made  to  operate  in  the  backfire  mode. 

The  backfire  bifilar  helical  antenna,  model  BBH-4,  has  a  diameter  of  1.77 
inches.  The  patterns  obtained  for  this  model,  when  a  series  of  brass  rods 
were  placed  inside  the  helix  and  coaxial  with  it,  are  shown  in  Figure  35, 

The  core  did  not  extend  beyond  the  feed  point.  The  core  diameters  ranged 
from  1/2  inch  to  1  1/2  inch  in  diameter  in  1/4-inch  steps.  The  frequency,  at 

which  these  patterns  were  taken,  is  one  hundred  megacycles  below  the  frequency 

■  '  ..  -  v>  .  '' 

given  bj^")£,(1.5  Gc)  for  this  model.  When  a  1  1/2-inch  core  is  used,  the  thick¬ 
ness  of  the  dielectric  polystyrene  foam  between  the  core  and  the  wires  is 
0,11  inches.  This  is  0.013  wavelengths  at  the  frequency  of  1,5  Gc. 
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Figure  35.  Patterns  of  BBH-4  With  a  Conducting  cbaxial  Core '  ^ 

at  f  =  1,4  G  ..  -  II 
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The  main  beam  of  the  patterns  shows  little  change  with  increasing  core 
diameter  except  for  some  asymmetry  due" to  misalignment  of  the  core  axis  and 

some  improvement  in  the  circular  polarization  on  axis.  The  principal  effect 

■  "  ..  .  '  .  -  •  ;  ..  "  .  >•  'r  .  ^  1  : 

of’ the  core  is  to  increase  the  excitation  of  the  higher  ord®£  Waveguide  mode. 

,  li  '\v  .? 

.This  is  seen  by  the  increase  ii5  endfire'  radiatiop.  The  increase  in  endfire 
radiation  is  significant  only  when  the  core  nearly  fills  the  center  of  the 
helix.  It  is  expected  that  a  further  increase  in  core  diameter  would  rapidly 

change  the  character  of  the  radiation  pattern,  but,  for  the  core  diameters 

.  a  "  ■■■  ”  ,( 

used  in  this  investigation,  surprisingly  little  effect  was  noted. 

Backfire  bifilar  helical  antennas  with  a  large  pitch  angle  tend  to  .  „ 

give  the  greatest  directivity.  To  provide  more  information  on  this  property, 
two  antennas  were  built  with  a  pitch  angle  of  45  .  BBH-8a  was  constructed  of 
number  twenty-two  wire  wound  on  a  1-inch  diameter  plexiglass  tube  48  inches 
long,  BBH-8b  was  constructed  of  number  thirty-four  wire  wound  on  a  0.475-inch 
diameter  pyrex  glass  tube  48  inches  long.  The  patterns  obtained  for  these 
antennas  are  shown  in  Figure  36. 

The  pattern  for  BBH-8a  at  1.315  Gc  has  a  28°  beamwidth.  The  length  of  the 

antenna  at  this  frequency  is  5.35  wavelengths.  Using  the  approximate  relation 

20 

between  beamwidth  and  gain  given  by  Kraus  ,  the  gain  of  the  antenna  at  this 
frequency  is  52.5  (17.2db).  This  gives  a  ratio  of  gain  to  length  in  wavelengths 

II 

of  9.82.  The  pattern  for  BBH-8b  at  2.575  Gc  is  20°,  and  the  length  of  the  antenna 

at  this  frequency  is  10.45  wavelengths.  Tjhe  gain  of  the  antenna  is  103.1 

(21.1db),  and  the  gain-to-length  ratio  is  9.9.  This  result  can  be  compared  with 

21 

the  expression  given  by  Schelkunoff  and  Friis  for  the  gain  of  a  long  array  of 

t»  ■  ’  '> 

isotropic  radiators  in  terms  of  its  length.. 

g  =  7.2  L/V 


’\\ 
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This  gain  results  when  Hansen-Woodyard  excess  phase  shift  is  used-in  the  array. 

^  -  If  '  ’  • 

The  backfire  bifllar  helix  provided  a  gain  about  1.375"  times  e's  great  without 
::  ...  .  „  -  ~ 
side  lobes.  If  side  lobes  were  permitted  as  shown"  for  ;))t he  1 owe#  frequencies,  ’ 

the  gain  is  increased  1.7  times.  The  backfire  bifilar  helices  are  .obviously" 

supergain  antennas  when  operated  near  the  critical  frequency.  The  usual 

instability  associated  with  supergain  afitennas  takes  the  form  of  a  very  rapid 

widening  of  the  beam  with  increasing  frequency.  One  half  percent  increase  in 

a 

frequency  increases  the  beam  width  2  . 

In  an  attempt  to  retain  some  part  of  the  high-gain  performance  of  the 
antenna  described  above  and  at  the  same  time  increase  the  band  width,  a  8-foot 
long  antenna  was  built,  tapered  from  one-half  inch  diameter  at  one  end  to 

■  t.  “  . 

one-inch  diameter  at^ the  other  end.  A  Rexolite  number  1422  core  was  used  with 
number  thirty- two  wire  wound  at  a  pitch  angle  of  45°.  The  patterns  obtained 
for  this  antenna  are  shown  in  Figure  37.  This  antenna  has  an  average  gain  of 
36(15. 5db)  in  a  five  percent  ban$  width.  Thus  tapering  the  helix  slows  the 
rate  with  which  its  patterns,  vary  with  frequency  and  at  the  same  time  reduces 

H 

the  maximum  gain  that  can  be  achieved.  This  is  the  result  expected  since  this 
antenna  corresponds  to  a  conical  log-spiral  antenna  with  a  very  small  cone  angle. 

A  backfire  monofilar  .helix  was  constructed ' to  show  that  the  nomofilar 
helix  will  also  operate  in  the  backfire  mode.  This  model  used  the  parameters  .. 
of  BBH-2.  The  conductor  in  this  case  was  a  Microdot  cable,  which  was  also 
used  to  carry  energy  from  the  feedpoint  to  the  receiver.  At  the  feedpoint, 
the  Microdot  cable  was  brought  radially  from  the  parameter  of  the  helix  to 
the  axis.  The  center  conductor  of  the  cable  was  connected  to  the  center  of  a 
straight  wire  ten-centimeters  long  placed  on  the  axis  of  the  helix,  The  patterns 
obtained  for  the  monofilar  helix  are  shown  in  Figure  38.  These  patterns  are  very 
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Figiire  37.  Measured  Patterns  for  a  Long  Tapered  Bifilar' Helix 


similar  to  those  obtained  for  the  bifilar  ^ptenna,  BBH-2.  This  demonstrates 

that  the  monofilar  helix,  will  operate  in  the  backfir'd  mode,  although  it  is 
./  ■  « 
more  difficult  to  food  than  the  bifilar  helix. 
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:7.  CONCLUSION  " 

The  object  of  this  research  was  to  provide  ^a^m^f^Saticaf  model  of  the 

"  ■■  ■  .  -  " 

backfire  bifilar  helical, .antenna  that  could  be  used  to  predict' the  radiation 

characteristics ^of  this  antenna .v  The  model  described  hore  was  shown  to  fulfill 
this  objective.  The  !^£udy  went  beyond  this  objective  in  providing  a  detailed 
experimental  study  of  the  radiation  patterns  of  the  antenna. 

In  this  research  the  radiation  patterns  Were  obtained  using  the  known 
determinantal  equation  for  the  helical  waveguide.  This  equation  was  factorized 
by  the  Wiener-Hopf  technique  ip  provide  the  Fourier  transform  of  the  current 
distribution  on  the  semi-infinity  helix.  The  radiation  pattern  was  related  to 
this  Fourier  spectruiy  *  A  /' 

Two  regions  in  the  frequency  spectrum  of  the  backfire  bifiiar  helix  were 

»  .  ■  » 

shown  to  be  important.  It  was  shown  that  at  the  lower  frequencies  the  bifilar 
helix  operates  as  a  waveguide.  Above  the  cutoff  frequency  of  the  principal 
waveguide  mode,  the  bifiiar  helix  operates  as  a  circularly  polarized  backfire 
antenna.  ' 

The  radiation  characteristics  of  this  antenna  were  shown  to  be  determined 
by  the  feed  region  current.  This  current  decays  with  distance  along  the 
structure.  The  near  field  has  a  phase  progression  directed  toward  the  feedpoint 
The  distance  to  which  the  feed  region  current  penetrates  the  structure  was  shown 
to  decrease  with  frequency. 

The  phase  pattern  of  the  antenna  was  not  obtained  in  the  calculations 
because  of  the " long  numerical  integration  involved  in  its  computation.  This 
prevented  a  more  accurate  accounting  of  the  endfire  radiation  from  the  antenna. 
The  good  agreement  between  measured  and  computed  results  reaffirmed  the  validity 
of  the  approximations  used  in  the  beginning  of  the  analysis. 


This  study  of  the  backfire  bifilar  helical  antenna  is  important  for 
several  reasons.  Not  lithe  least  of  these  is  the  fact  that  the  antenna  is  a 

useful,  practical  device  in  itself.  It  provides  circularly  polarized  fields 

/  ;  "  ■ '  .. 

with  a  single  lobed  pattern  of  almost  any  desired  gain,  When  iiitch  angles 

...  ■ '  "  "  •/ 

near  45°  are  used,  the  antenna  can  provide  circularly  polarized  fields  with  a 
"!■■■'  ■  ...  '  v 
pattern  that  is  omnidirectional  "in  azimuth  with  high  vertical  directivity. 

If  i^/is  inconvenient  to  provide  a  balanced  feed  for  this  antenna,  it  can  be 

fed  from  a  Coaxial  line  run  inside  one  of  the  helical  conductors  to  the 

feedpoint  on  the  axis.  The  decay  of  the  feed  region  current  will  keep  antenna 

■  currents  from  flowing  on  the  unbalanced  feed  line. 

This  study  provides  further  insight  into  the  operation  of  the  conical 

log-spiral  antenna,  and  it  is  the  first"  simply  periodic  backward  wave  or 

backfire  antenna  for  which  a  theoretical  solution  has  been  obtained.  The 

.  "  .  -  c, 

study  of  other  types  of  periodic  antenna  structures  that  operate  in  the 
backfire  mode  will  aid  in  understanding  the  operation  of  known  log-periodic 
antennas  and  will  suggest  new  ones.  These  structures  are  important  in 
themselves  for  their  unidirectional  single  lobe  patterns  and  their  frequency 
scanning  characteristics,  as  well  as  for  the  insight  they  give  in  the  operation 
of  log-periodic  antennas. 

"  The  existence  of  a  decaying,  backward  wave  current  in  the  feed  region 
suggests  that  the  patterns  of  the  backfire  helix  can  be  explained  in  terms 
of  improper  (decaying)  helical  waveguide  modes.  An  investigation  along  these 
lines  is  now  in  progress  by  another  investigator.  The  results  of  such  an 
analysis  may  be  simpler  in  interpretation  than  the  present  analysis. 

a  -  *  ••  w 

Another  line  of  inquiry  suggested  by  this  work  concerns  the  effect  of  a 
conducting  core  on  the  axis  of  the  helical  antenna.  In  this  study  it  was  shown 


that  the  helical  conductors  could  be  placed  remarkably  close  to  the  conducting 
core  with  little  degradation  of  performance.  This  study  should  be  extended  to 

X  .  •'  C  ■  • 

other  helix  parameters.  The  effect  of  extending  the  conducting  Core  beyond  the 
feedpoint  should  be  investigated.  A  backfire  monofilar  helix  fed  against  a 
conducting  central  core  should  also  be  investigated. 

/(The  study  of  long  slowly  tapered  bifilar  heliceis  should  be  extended.  The 
results  of  this  Study  combined  with  results  obtained  for  the  conical  log-spiral 
antenna  would  establish  Jhe  relation  between  d'ain,  band  width,  and  size  over 
the  entire  spectrum  of  helical  or  conical  log-spiral  parameters.  The  conical 
log-spiral  antenna  is  known  to  have  practically  constant  characteristics  over  %• 
ten-to-one  or  greater  band  widths.  This  antenna  in  the  .usual  range  of  parameters 

'  ...  .  .  ‘  I,  / 

has  considerably  less  gain  than  can  be  obtained  by  a  slowly  tapered  helix,  while 

'  ■  ..  '  /  ■ 

the  slowly  tapered  helix  has  a  relatively  smaller  band?width, 

It  is  characteristic  of  a  study  of  a  new  device  that  many  diverging  courses 

("■■■■ 

of  investigation  are  suggested.  A  few  of  these,  that  were  suggested  during  the 

'  '  !l 

study  of  the  backfire  bifilar  helical  antenna,  are  described  above.  However, 
the  major  effort  in  this  investigation  was  concentrated  on  the  study  of  the 
radiation  pattern- of  the  antenna.  In  this  w^y  it  has  been  possible  to  provide 
information  on  this  aspect  of  the  problem  over  on  extensive  range  of  parameters. 
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APPENDIX  A 


A  FOURIER  TRANSFORM 


The  Fourier  transform  of  the  free  space  Green's  function  specialized  to 
the  helix  is  important  to  the  development  of  the  analysis  of  the  backfire 
bifilar  helical  antenna.  This  transform  is  easily  related  to  the  transform 


of  the  function 


<z).  = 


o  2  2  1/2 

exp(jk[A  -2AB  cos  (c|>tTz)  ■»  B  +  (z-d)  ]  _  ) 

1/2 


2  2 
[A  -2AB  cos  (4h-tz)  +  B  +  (z-d)] 


This  function  is  periodic  in  u  =  <j>fTz  and  carl  be  expanded  in.  the  Fourier  series 


1  y  -  jn(4>+Tz) 

gl(z)  =2 i  E  e 

n=-oo 


2  22  1/2 
exp(jk[A  -2AB  cos  u-t-B  +(z-d)  ]  )  eJnudu  (A_2) 

[A2-2AB  cos  u+B2+ (z-d) 2 ] 


If  we  take  the  Fourier  transform  (Equation  A-2)  and  change  the  order  of  inte¬ 
gration,  the  resulting  expression  is 


~  1  y  "jn<]> 

g1(B)=  — g  ^  e 

4ir  n=-oo 


1J"£ 


exp(jk[A2-2AB  cos  u+B2+(z-d)2] 1//2)  j(B-nT)z,  , 
- : - - -  e  dzdy 

2  2  o  * 

[A  -2AB  cos  u+B  +(z-d)  ] 


Using  the  result  listed  in  Magnus  and  Oberhi'ettinger  22  this  expression  becomes 


g.(3)  = £  e'Jn<})  /  eJnUH(1)  (rk2-(B-nT)2]1/2[A2-2AB  cos  u+B2]1/2ej(8"nT)d. 

1  oil  I  O  uV 

n=-oo  Jo  J 

The  addition  theorem  tor  the  Hankel  function,  applied  to  the  expression,  gives 


■l.«)-  fc  L  9 


27T 

-*H/£  ; 

I  n--K 

Jo 


j  [(B-n)d-n<f>]  f  L  j  (B|k2_(B_njT)2]V2)Ha)(A[k2_(3_mT)2]V2)ej(m+n)udu 
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Interchanging  the  order  of  integration  with  summation  and  using  the  orthogonality 

*.\  ••  ...  ..  ,  ..  ,Ji  • 

of  the  exponential  functioii'  allows  the  result  to  be  written 


gx(B)=  J'E  Jm(6[k2-(B-mTlj2]1/2)H^1U['k2-(6-mT)2]1/2)  ed  (S+mT^0Jm ^ 


This  expression  is  valid  wher  B  is  less  than  A.  If  A  is  smaller  it  is 

necessary  to  interchange  A  and  B.  If  A  equals  B  the  series  does  not  converge. 

23  \\  •  « 

A  similar  result  was  obtained  by  Phillips  another^way .  Hejj  gives  an  extensive 
discussion  of  the  convergence  of  this  series.  | 


If  g^  is  deduced  from  in^Equation  (Arl)  by  changing  the  sign  of  the 

i.  "  "  /v 

cosine  term,  g^  can.  be  deduced  from  g^  by  changing  the  sign  of  all  odd  numbered 


g  (13)  »  J  E  <-l)m  J  (B|ik2-(B-mT)2]1/2)  H(i)  (A[k2-(B-nr)2]1/2) 
-5^  m  m 


j(B+mT)d  jm<j> 
e  e 


These  transforms,  then,  have  the  sum 


g  (B)+g  (B)  =  §E  J  (B^2-(B-mT)2]1/2)  H(1)(A[k2-(B2-mT)2]1/2) 
&  ±  ^  m  m 


i(B+mi)d 


and  the  difference 


g  (B)-g  (8)  E  J  (B^2-(6-PT)2]1/Z)  HCl)(A[k2-(B-pT)2]1/2) 

±  £t  £$  P  P 


p=2ti+l 


„j(lfc-pT)d  jptf’ 


r 


-4 


"ft 
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)/  Finally,  the  Fourier  transform  of  the  complex  conjugate  of  g(z)  is  deduced  from 

•  is*  .  ;;  "  -  Cl) 

g(3)  by  replacing  the  Hankel  function  of  the  first  kind,  H  with  the  Hankel 

X  '  ..  '  n-  - 

(2) 

function  of  the  second  kihd.H  . 

’  n  \»  ■ 


no 


APPENDIX  B 


THE  RELATION  BETWEEN  THE  FOURIER  SPECTRUM  OF  A  LINEAR 
CURRENT  DISTRIBUTION  AND  ITS  RADIATION  PATTERN 

The  vector  potential  of  a  linear  current  distribution  along  the  z-axis 


is  giVeh  by 


A(r)  =  u  I  g(r,S;.z)  I  ^  dz 


A 


where  u  is  the  direction  of  the  current  elements  and 


CB-1) 


g  = 


exp(-jk[(r  sin  6)2+(z-r  cos  9)2]^2) 

2  2 
47f[(r  sin  9)  +(z-r  cos  9)  ] 


By  Parseval’s  theorem  this  is  equivalent  to 


A (r)  =  u  21T  l  gT  (r,9;|J)  f(B)  dB 


(B-2) 


where 


KB)  =  I  Kz)  ejBz  dz. 


and 


G(r,  9;  13)  =  ^ 


g*(r, 9; z)e^Bzdx 


This  is  the  same  Fourier  transform  pair  used  in  Appendix  A.  In  this  case. 


j  ..i-d)  .  Or,2  n2l1/2\  C0S  6 

G(r,  ;  13)  =  ~  H  (r  sin  Ok  -B  ]  )  e 

’  Sir  o 


Thus  Equation  (B-2j  is 
00 


A(r)  = 


-ju 


/ 

-00 


H^2) (r  sin  9[k2-B2]1/2)  1(B)  e“jSr  COS  GdS 


(B-3) 


To  obtain  the  radiation  pattern  from  this  expression,  it  is  necessary  only  to 
obtain  an  asymptotic  estimation  f.or  the  integral  yalid  for  large  r.  This  is 


f's 


i 
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We  have  deduced  by  a  long  and  complicated  argument  the  well  known  result  that 

•  ••  *  ’’ 

the  current  distribution  of  a  line  source  and  its  radiation  pattern  are  Fourier 

transform  pairs.  This  result  is  more  easily  obtained  by  assuming  that  r  is  much 
••larger  .than..* -in  .Equation  (B-l)  giving  „  " 


■  J 

j' 


A(r)  ~ 


“-jkr 
e  J 

4ffr 


os 

/ 


jk  cos  S  i 
e  .I(z)  d 


This  more  direct  approach,  however,  is  valid  only  when  the  current  has  finite 
support.  If  the  current.  hasAfinit.e  support  and  is  square  integrable  there, 
•Theorem  X  of  Paley  and  Wiener*  stages  that  its  Fourier  spectrum  is  entire. 

The  current  spectrum  then  will  contribute  no  singularities  of  its  own  tC  P(u) 
in  Equation  (B-6) .  If  the  current  distribution  is  not  finite,  its  spectrum 
may  have  singularities  of  its  own  that  must  be  considered  in  the  evaluation  of.. 
Equation  (B-6) . 

Mathematically,  such  additional  singularities  can  be  handled  by  the  method 
of  Van  der  Waerden'1''1' .  Physically,  they  are  usually  interpreted  as  surface 
waves  guided  by  the  structure.  The  contribution  of  the  branch  point  at  u  =  j 
might  be  called  the  "spa^e  wave"  of  the  distribution.  It  is  this  "square  wave" 
that  we  have  identified  with  the  radiation  patterns  of  a  current  distribution 
in  this  study. 
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